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The evolution of the Universe between inflation and the onset of Big Bang Nucleosynthesis is 
difficult to probe and largely unconstrained. This ignorance profoundly limits our understanding of 
dark matter: we cannot calculate its thermal relic abundance without knowing when the Universe 
became radiation dominated. Fortunately, small-scale density perturbations provide a probe of 
the early Universe that could break this degeneracy. If dark matter is a thermal relic, density 
perturbations that enter the horizon during an early matter-dominated era grow linearly with the 
scale factor prior to reheating. The resulting abundance of substructure boosts the annihilation rate 
by several orders of magnitude, which can compensate for the smaller annihilation cross sections 
that are required to generate the observed dark matter density in these scenarios. In particular, 
thermal relics with masses less than a TeV that thermally and kinetically decouple prior to reheating 
may already be ruled out by Fermi-LAT observations of dwarf spheroidal galaxies. Although these 
constraints are subject to uncertainties regarding the internal structure of the microhalos that form 
from the enhanced perturbations, they open up the possibility of using gamma-ray observations to 
learn about the reheating of the Universe. 


I. INTRODUCTION 

The expansion history of the Universe prior to the on¬ 
set of Big Bang Nucleosynthesis (BBN) is unknown, but 
the nearly scale-invariant spectrum of the primordial cur¬ 
vature perturbations provides strong evidence that the 
Universe experienced a period of inflation EH shortly 
after the Big Bang. Inflation leaves the Universe devoid 
of radiation, and yet the primordial abundance of light 
elements indicates that the Universe was radiation dom¬ 
inated during BBN. Therefore, the Universe must have 
become radiation dominated at some point after infla¬ 
tion, but we do not know when that transition occurred, 
nor do we know how the Universe evolved between infla¬ 
tion and the onset of radiation domination. The existence 
of the cosmic neutrino background provides the most ro¬ 
bust constraint on the temperature of the Universe when 
it became radiation dominated: this reheat temperature 
must exceed 3 MeV to generate the neutrinos required 
to produce the observed abundances of light elements 
[3-0| and the observed power spectra of anisotropies in 
the cosmic microwave background and large-scale den¬ 
sity perturbations The vast difference between 3 

MeV and the theorized energy scale of inflation leaves 
a tremendous gap in our understanding of the thermal 
history of the Universe. 

In the simplest models, inflation is powered by a single 
scalar field, the inflaton. Inflation ends when the inflaton 
begins to oscillate around the minimum of its potential, 
and the Universe becomes radiation dominated when the 
inflaton eventually decays into relativistic particles E3- 
E3- If the potential is quadratic, a coherently oscillat¬ 
ing scalar field has the same dynamics as a pressureless 
fluid EH- Therefore, most inflationary scenarios include 
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a transient period of effective matter domination between 
the end of inflation and radiation domination. Further¬ 
more, after the inflaton decays into radiation, other os¬ 
cillating scalar fields may come to dominate the energy 
density of the Universe; such scalars commonly result 
from stabilized moduli in string theories EB“I13 
responsible for generating the primordial curvature per¬ 
turbation in the curvaton model [23 - [^ . It is also pos¬ 
sible that quasi-stable massive particles could dominate 
the Universe prior to BBN [e.g. [l^. Thus, a pre-BBN 
early-matter-dominated era (EMDE), driven by either an 
oscillating scalar field or a quasi-stable massive particle, 
is a generic prediction of several early-Universe theories. 

There may have been several EMDEs, as the energy 
density of longer-lived oscillating scalar fields could sur¬ 
pass the energy density of the relativistic decay prod¬ 
ucts generated at the end of earlier EMDEs. In this 
case, the last EMDE is most cosmologically significant 
because it generates the current content of the Universe 
and its end determines the timing of reheating, which is 
defined to be the onset of the final radiation-dominated 
era. Reheating is usually assumed to occur too early 
to have any cosmological impact, but there are strong 
motivations to consider lower reheat temperatures. The 
moduli fields predicted by string theory have gravita¬ 
tional couplings, so their masses must be .^10 TeV to 
avoid violating the BBN bound on the reheat tempera¬ 
ture (Trh ^ 3 MeV). String theories also generally pre¬ 
dict that the lightest modulus field should not be signif¬ 
icantly more massive than the gravitino, which implies 
that the modulus mass should be less than ^1000 TeV if 
supersymmetry is to mitigate the electroweak hierarchy 
problem. Consequently, the reheat temperature in these 
theories is generally less than a few hundred GeV (see 
Ref. [1^ for a recent review). 

Our ignorance of the pre-BBN thermal history pro¬ 
foundly limits our understanding of the origins of dark 
matter [e.g. the origins of the baryon asymme- 
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try [e.g-HO, m, and the connection between the pri¬ 
mordial power spectrum and the inflaton potential [e.g. 
Idrl - f^ . In particular, the relic abundance of dark matter 
is profoundly altered if dark matter thermally decouples 
during an EMDE [29l - l3li the decaying field di¬ 

lutes the relic abundance of thermal dark matter, so a 
smaller value of {av) is required to generate the observed 
dark matter density. The dark matter density may be 
enhanced, however, if the dark matter is also produced 
nonthermally (as a decay product, for instance), in which 
case {av) may need to be increased so that the excess 
dark matter is eliminated. Ref. [d^ demonstrated that 
nearly any supersymmetric dark matter particle, with a 
wide range of {av) values, can give the correct relic abun¬ 
dance for some combination of a low reheat temperature 
and nonthermal production. Even in the absence of non- 
thermal production, it is possible to obtain the observed 
density of dark matter if {av) <C 3 x 10“^® cm^ s“^. 
Therefore, constraints on the dark matter annihilation 
rate cannot rule out a thermal origin for dark matter as 
long as the reheat temperature remains unconstrained. 


In this work, I show that an EMDE’s effect on density 
perturbations on scales that enter the cosmological hori¬ 
zon prior to reheating can significantly enhance the dark 
matter annihilation rate by increasing the abundance of 
microhalos. During an EMDE, subhorizon perturbations 
in the dominant energy density grow linearly with the 
scale factor [dj |d^ . If the duration of the EMDE is suf¬ 
ficiently long, these perturbations become nonlinear and 
could produce gravitational waves [d^, , but probably 

not black holes [2^. Ref. [5l| demonstrated that non- 
relativistic decay products produced during the EMDE 
inherit these enhanced perturbations. Consequently, if 
dark matter is generated nonthermally and is nonrela- 
tivistic at reheating, perturbations in the dark matter 
density that enter the cosmological horizon during the 
EMDE are significantly enhanced compared to larger- 
scale perturbations. Ref. showed that this small- 
scale inhomogeneity would radically increase the number 
density of sub-Earth-mass microhalos. Unfortunately, 
these microhalos are too small to detect with gravita¬ 
tional lensing in the near future (s^ . This analysis was 
later extended to scenarios in which the dark matter is 
generated nonthermally and then annihilates down to the 
observed abundance [53, k 


___ j|. In these scenarios, sub¬ 
horizon perturbations in the dark matter density grow 
linearly during the EMDE and then decrease slightly at 
reheating. Although this brief decrease does not erase 
the growth of the perturbations during the EMDE, it is 
likely that density perturbations on these small scales are 
further suppressed by free-streaming of the dark matter 
particle after reheating M|. 

However, if the dark matter is produced thermally 
during an EMDE, it may kinetically decouple prior to 
reheating, in which case the free-streaming scale is far 
smaller than the horizon scale at reheating ^ . Further¬ 
more, the dark matter’s thermal origins demand that it 
must annihilate in high-density environments, which pro¬ 


vides a means of detecting the microhalos generated by 
the EDME. This scenario is the focus of this work; I 
will assume that the dominant energy component dur¬ 
ing the EMDE either does not decay into dark matter 
or that the branching ratio for decay into dark matter 
is small enough that the generated dark matter density 
is far smaller than the thermal relic density. I will show 
that subhorizon perturbations in the thermally generated 
dark matter density grow linearly during EMDE and that 
the resulting matter power spectrum is identical to the 
power spectrum derived in Ref. (Fij i for nonthermal dark 
matter in the absence of annihilations. 

If the dark matter kinetically decouples prior to re¬ 
heating, the largest scales enhanced by the EMDE are 
not suppressed by the free-streaming of the dark mat¬ 
ter particle nor by its elastic scatterings with relativistic 
particles. Therefore, these enhanced small-scale pertur¬ 
bations will significantly increase the abundance of mi¬ 
crohalos: instead of being rare, isolated objects, micro¬ 
halos contain most of the dark matter in these scenarios. 
Even though these particles must have annihilation cross 
sections that are far smaller than the canonical thermal 
cross section, the boost to the annihilation rate from the 
microhalos is sufficient in some cases to increase the an¬ 
nihilation rate beyond current observational limits from 
the Fermi Gamma-Ray Telescope [FFl [F^- This result 
demonstrates that gamma-ray observations are capable 
of probing the evolution of the Universe prior to BBN 
and constraining the origins of dark matter. 

I begin in Section IB Al by reviewing the thermal pro¬ 
duction of dark matter during an EMDE. In Section lllBI 
I analyze the evolution of density perturbations during 
the EMDE and after reheatin g, r ecovering the matter 
power spectrum found in Ref. The suppression of 

small-scale perturbations due to elastic scattering and 
free-streaming of the dark matter particles is considered 
in Section uni This section also explores the impact the 
decoupling scale has on the microhalo population. In Sec- 
tion llV Al I estimate how these microhalos boost the dark 
matter annihilation rate, and I compare the resulting an¬ 
nihilation rate to constraints from gamma-ray observa¬ 
tions in Section IIVBI Finally, I summarize my results 
and discuss avenues for future investigation in Section 
El The appendices present the derivation of the pertur¬ 
bation evolution equations and their initial conditions. 
Natural units {h = c = ks = ^) are used throughout this 
work. 


II. THERMAL DARK MATTER DURING AN 
EARLY MATTER-DOMINATED ERA 

A. Background Evolution 


The dominant component of the Universe prior to re¬ 
heating is modeled as a pressureless fluid that decays 
perturbatively into relativistic particles. Although I will 
refer to this pressureless fluid as a scalar field, it could 
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equally well be a collection of massive particles. I as¬ 
sume that dark matter is created thermally and that no 
dark matter is produced during the decay of the scalar 
field. In this scenario, the background equations for the 
energy density of the scalar field, the radiation den¬ 
sity pr (not including dark matter particles, even if they 
are relativistic), and the number density of dark matter 
particles are 

"k pif, = (1^) 

^Pr + AHpr = T^p^ + {crv){E^) {n\ - , (lb) 

+ =-{av) , (Ic) 

where is the scalar decay rate; {av) is the velocity- 
averaged cross section for dark matter annihilations; 
{E-^ = Pxl^x is the average energy of a dark matter 
particle; and is the number density of dark matter 
particles in thermal equilibrium.^ For a dark matter par¬ 
ticle with mass and internal degrees of freedom. 
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where T is the temperature of the radiation bath. When 
the dark matter particle is nonrelativistic, {E^) ~ m-)^, 
while {E^) ~ 3.15T if rrix <C T. To facilitate a numerical 
solution to the background equations, I make the ap¬ 
proximation that {E^) ~ -I- 9.93T^, which matches 

Px/n^ to within 10%. This system of equations was stud¬ 
ied extensively in Ref. [30l |: in this section I will only 
summarize the key results that are required to compute 
the relic abundance of dark matter in these scenarios. 

Figure [T] shows the evolution of the densities deter¬ 
mined by solving Eq. ©• The scalar field dominates the 
energy content of the Universe until F^ ~ iF. At that 
time, p 0 begins to decay exponentially, and the Universe 
becomes radiation dominated. Although radiation does 
not instantaneously dominate the energy density when 
F^ = iF, it is still useful to define the reheat tempera¬ 
ture Trh by the relation 
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where mpi = is the Planck mass, and 

g*{T) = pr.(T)/[(7r^/30)T^] is the number of relativistic 


^ These equations differ slightly from those in Ref. [^ . which in¬ 
cludes an extra factor of two in the radiation equation to account 
for the fact that two photons are produced during each annihila¬ 
tion event. This factor of two is not necessary if the dark matter 
is composed of Majorana particles, because then the rate of an¬ 
nihilation events per volume is ((Tu)n^/2. However, its inclusion 
does not affect the relic abundance of dark matter. 



FIG. 1: The densities of the decaying scalar field, radiation, 
and matter during reheating. In this hgure, nix = 750 GeV, 
and H = Ttj, when ajai — “2 y, 10®, which corresponds to a 
reheat temperature of 5 GeV. The two solid curves show px for 
the two values of {av) that generate the dark matter density 
observed today {Q.xh? = 0 . 12 ) for these values of and Trh. 
The top solid curve shows a freeze-out scenario with {av) = 
9.0 X 10“®^cm®s“^, whereas the bottom solid curve shows 
a freeze-in scenario with {av) = 2.9 x 10“®®cm®s~^. The 
dotted curve shows the equilibrium density of dark matter, 
Peq = {Ex)nx,eq- 


degrees of freedom. I evaluate g*{T) by summing the 
contributions of all Standard Model particles to Pr (see 
Appendix A of Ref. [H^). This definition of the reheat 
temperature can be used to eliminate F^ from Eq. O; 
Trh then parametrizes the onset of radiation domina¬ 
tion. It is also useful to define urh to be the value of the 
scale factor when F^ = iF. Since El oc during the 

EMDE, 
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where a/ is the value of the scale factor at some point 
during scalar domination, p^j is p^ when a = o/, and I 
have introduced the shorthand 5 *rh = g*(TRH)- In Fig- 
ure[Tl Trh = 5 GeV and aRu/a/ = 2 x 10®. 

Prior to the onset of radiation domination, the trans¬ 
fer of energy from the scalar field to radiation forces 
Pr oc This scaling can be derived from Eq. m 

in the limit that T^p^ dominates over all other contribu¬ 
tions to dpr/dt] since oc and iF oc during 

the EMDE, Eq. (fib)) then implies that p(.(a) oc To 

obtain an explicit expression for pr during the EMDE, 
however, the AHpr term in Eq. m must be included, 
and an initial value of pr must be provided. If pr = 0 
when a = aj, then 


Pr — 


2^.2 / ii^5*RhP0,/ 

5 30 



- 3/2 
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As shown in Figure [1] pr increases rapidly after a = aj, 
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FIG. 2: The present-day abundance of dark matter as a function of the velocity-averaged annihilation cross section (av). 

For small (av) values, the dark matter “freezes in” and 12^ oc (uv), whereas for larger (av) values, the dark matter “freezes 
out” and oc (av)~^. The left panel shows that increasing decreases in both cases. The right panel shows 

that decreasing Trh while keeping m^/TaH fixed also decreases albeit mildly if (av) is in the freeze-out regime. In both 

panels, the curves are based on numerical solutions to Eq. (Hj, while the open circles show the predictions of Eqs. © and (f8)l. 
The dashed line marks the observed value established by the Planck mission: — 0.12 [s^ . 


reaches a maximum when a/ai ~ 1.48, and then de¬ 
creases as pr oc until the onset of radiation dom¬ 

ination. Since the radiation cools during the EMDE 
(T cx Trh is not the maximum temperature of 

the radiation bath, and it is possible to thermally pro¬ 
duce particles with ^ Trh. Dark matter is efficiently 
produced provided that < Tmax, where 

Tmax „ [mpig(aj)]V'^ 

Trh - ■ g*(TMAx)i/4 ^ ’ 


As long as Tmax 2mx, the relic abundance of dark 
matter will be independent of Tmax, which implies that it 
is also independent of H(ai). Throughout this analysis, 
I will chose p^j such that Tmax = Sm^, which ensures 
that the results are not sensitive to the details of the 
onset of the scalar decay. 

There are two thermal dark matter production sce¬ 
narios illustrated in Figure [T] that both produce the 
same present-day dark matter density. If (av) is suffi¬ 
ciently large that pair production brings the dark mat¬ 
ter into thermal equilibrium, then = n^,eq while 
H < {av)n^^eq- When H ~ {av)n^^eq, the dark matter 
“freezes out,” and cx a~^ thereafter. If dark matter 
freezes out prior to reheating. 


~ 1.6 X 10 
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where Tf is the freeze-out temperature defined by 
H(Tf) = {av)n^^eq{Tf). The ratio m^/Tf depends only 
logarithmically on (av), Trh, rnost cases 

5 ^ m^/Tf ^ 25. Equation © indicates that increas¬ 
ing TOx/Trh dramatically reduces the relic abundance of 


dark matter, as illustrated in Figure [H Since m^/T/ is 
nearly constant, increasing TOx/Trh is roughly equivalent 
to decreasing Trh/T/. The resulting decrease in 
has a simple explanation: if the dark matter freezes out 
during the EMDE, then the subsequent transfer of energy 
from the scalar to radiation will increase the number of 
photons per dark matter particle, which dilutes the relic 
abundance of dark matter. As shown in Figure [H the ob¬ 
served abundance of dark matter = 0.1197±0.0022 

[ 5^ 1 can only be generated if {av) <C 3 x 10“^® cm® s“®. 

As {av) decreases, the relic abundance of dark mat¬ 
ter increases, until {av) becomes small enough that pair 
production cannot bring up to its equilibrium value 
before pair production ceases when T ^ my-jA. At that 
point, the dark matter “freezes in,” and the comoving 
number density of dark matter is subsequently preserved, 
as demonstrated in Figure [T] If the dark matter becomes 
nonrelativistic prior to reheating, the freeze-in abundance 
of dark matter is 


~ 0.062 


3/2 

g*RH 
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In contrast to the freeze-out scenario, where increasing 
{av) lowers the dark matter abundance by making dark 
matter annihilation more efficient, annihilations are not 
common in the freeze-in scenario, and increasing {av) 
increases the relic abundance by making pair produc¬ 
tion more efficient. However, increasing mx/TRH still 
decreases the relic abundance in the freeze-in regime for 
the same reason as before: as TOx/Trh increases, more 
photons are created after the production of dark matter 
halts. 
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The final abundance of dark matter in both the freeze- 
in and freeze-out scenarios was derived in Ref. [s^ using 
the approximation that the transition from matter dom¬ 
ination to radiation domination occurs instantly when 
T = Trh- These analytic estimates accurately recover 
the functional dependence of the dark matter density on 
’Trh, and (crv), but their normalization must be ad¬ 
justed slightly to correct for the fact that reheating is 
not instantaneous. The expressions for given by 

Eqs. ([7|) and (|S]) have been normalized to match the dark- 
matter density obtained by numerically solving Eq. O- 
Figure [5] demonstrates how well these expressions agree 
with the numerical results for a wide range of parameters. 


B. Perturbation Evolution 

An EMDE also affects the evolution of the density per¬ 
turbations in both the radiation and the dark matter 
[ 5 ^ . The evolution of perturbations in the mat¬ 

ter density is primarily determined by the expansion rate 
of the Universe. While the Universe is radiation domi¬ 
nated, the fractional perturbation in the matter density 
grows logarithmically with the scale factor after the per¬ 
turbation mode enters the horizon (fc > aH, where k is 
the wave number of the plane-wave perturbation). When 
the Universe becomes matter dominated, subhorizon per¬ 
turbations begin to grow linearly with the scale factor. 
An EMDE has the same effect on the dark matter density 
even though a different pressureless fluid dominates the 
Universe. If dark matter is generated during scalar de¬ 
cays and does not annihilate, density perturbations in the 
dark matter on subhorizon scales grow linearly with the 
scale factor during an EMDE, and then they smoothly 
transition to logarithmic growth when the Universe be¬ 
comes radiation dominated [5l|. 

Subhorizon density perturbations in the scalar field 
also grow linearly with the scale factor during the EMDE. 
The relativistic decay products of the scalar field inherit 
this inhomogeneity; subhorizon perturbations in the radi¬ 
ation density grow during an EMDE, but that growth is 
lost when the Universe becomes radiation dominated [sij . 
When dark matter annihilations are added to scenarios 
in which dark matter is primarily produced during scalar 
decays, the dark matter perturbations still grow linearly 
during the EMDE, but they decrease slightly when much 
of the dark matter annihilates during the transition to 
radiation domination . 

In this section I explore how an EMDE affects pertur¬ 
bations if the dark matter is not produced in scalar de¬ 
cays and thermally decouples prior to reheating. Evolu¬ 
tion equations for the fractional density perturbations in 
p^, p^, and pr [di = {p%- p\)lp\-, where p\(t) is the spa¬ 
tially uniform background density] are obtained by per¬ 
turbing the covariant form of the energy-transfer equa¬ 
tions given in Eq. ©■ This procedure and the resulting 
evolution equations are presented in Appendix The 
annihilation terms in Eq. induce a coupling between 



FIG. 3: The evolution of the density perturbation in dark 
matter for a mode that enters the horizon prior to reheat¬ 
ing {k = 35feRH,h) in the freeze-in scenario {{cv) = 2.9 x 
10”®^ cm® s~^). The background evolution for this scenario is 
shown in Figure[T]with aj = 3.6 x 10“"^. The solid curve shows 
Sx- The long-dashed curved shows the pre-reheating evolution 
of S^, and the short-dashed curve shows the post-reheating 
evolution of Sx predicted by Eq. uni) with aiog = 1.29aRH. 


Sx and Sr, and these interactions are included in the per¬ 
turbation equations. Elastic scatterings between dark 
matter particles and relativistic particles lead to addi¬ 
tional couplings between Sx and Sr', these effects are con¬ 
sidered in Section IIIII along with the suppression of per¬ 
turbations due to the free-streaming of the dark matter 
particles. 

The perturbation evolution equations are solved nu¬ 
merically for each k value starting before the mode enters 
the horizon and after the dark matter particle becomes 
nonrelativistic. Since the scalar field is the only source of 
energy at the start of the EMDE, the perturbations are 
adiabatic and are all linearly related to the initial per¬ 
turbation in the spacetime curvature $ 0 ^ see Appendix 
m for details. Figure [3] shows the evolution of the dark 
matter density perturbation Sx for a mode with wave 
number k = 35fcRH,hj where fcRH.h = aRHl^(aRH) is the 
wave number of the mode that enters the horizon at re¬ 
heating. The background evolution for this cosmological 
scenario is shown in Figure [TJ the value of o/ in this fig¬ 
ure is determined by setting a = 1 at the start of the 
numerical integration of the perturbation equations. In 
Figure[3l (av) = 2.9x 10“®® cm® s“^, and the dark matter 
freezes in; the comoving number density of dark matter 
particles is conserved for a ^ 10, which corresponds to 
a/ai ^ 3x 10"^ in Figure[TJ When pair production halts, 
the mode shown in Figure [3] has not yet entered the hori¬ 
zon. Since both p^ and Px are proportional to a“® at this 
point onwards, adiabaticity demands that Sx = S^ on su¬ 
perhorizon scales. The long-dashed line in Figure|3]shows 
the evolution of S^ prior to reheating [see Eq. (|Blbl) in 
Appendix IB] : we see that Sx S^ when a ~ 10. 

The evolution of the dark matter perturbation in the 
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FIG. 4: The evolution of the density perturbation in dark matter for modes that enter the horizon prior to reheating {k = SSIcRH.h 
in the left panel and k = 95fcRH,h in the right panel) in the freeze-out scenario ((crv) = 9.0 x 10“®^ crn^s”^). The background 
evolution for this scenario is shown in Figure[T]with a/ = 3.4 x 10“^. The solid curve shows 3^. The long-dashed curved shows 
the pre-reheating evolution of S^, and the short-dashed curve shows the post-reheating evolution of 3^ predicted by Eq. (HSl) 
with aiog = 1.29aRH. 


freeze-out scenario is shown in Figure S) While the 
dark matter is in thermal equilibrium, 6^ ~ i5^,eq, where 
i55(._eq = n'x,eq/'^x,eq ~ 1 IS the perturbation in the equi¬ 
librium density of the dark matter. Since n^^eq is de¬ 
termined by the radiation density, (^^.eq is related to the 
perturbation in the radiation density: 

<5x.eq = (9) 

The evolution of (5x,eq is also shown in Figure ID While 
the mode is outside the horizon, 6r is constant, but (5x,eq 
increases as m^/T increases. After the mode enters the 
horizon, Sr begins to grow as it is sourced by the increas¬ 
ingly inhomogeneous decaying scalar field. However, the 
amplitude of Sr plummets during reheating, and after 
reheating, Sr oscillates with a constant amplitude that 
is far smaller than it would have in the absence of an 
EMDE. This strange behavior of the radiation perturba¬ 
tion was discovered in Ref. (5lj | and later confirmed in 
Refs. Since d^.eq — 0.25(m^/r)i5r while the dark 

matter is nonrelativistic, d^.eq ^Iso grows prior to reheat¬ 
ing and then decreases rapidly before beginning oscilla¬ 
tions, as seen in Figure SI If the dark matter freezes out 
(or kinetically decouples) after reheating, this reduction 
in Sr would further suppress the small-scale dark matter 
density fluctuations. However, it is doubtful that this 
additional suppression on top of the usual suppression 
on scales smaller than the decoupling scale would be ob¬ 
servable. Instead, I will focus on scenarios in which the 
dark matter freezes out and kinetically decouples prior 
to reheating. 

Figure S] shows the evolution of two such perturba¬ 
tion modes, one that freezes out prior to horizon en¬ 
try (fc = 35fcRH,h) and one that freezes out after hori¬ 
zon entry {k = 95fcRH,h)- In both cases, S^ ~ '^x.eq 
until H = (avjn^^eq, which is marked as “Freeze-out” 


in Figure m After freeze-out, S^ matches the evolution 
of S^, as it did in the freeze-in case shown in Figure |31 
In the absence of pair production and annihilations, the 
scalar perturbations and the dark matter perturbations 
are governed by the same set of equations prior to re¬ 
heating: see Eq. (IA18I) . The solution to these equations 
is S(a) = So + (2/3)fc^<i>on, where do is a constant and 
k = k/H{a = 1). Therefore, after freeze-out or freeze-in 
and before reheating, S^ and S^ differ by a constant, and 
that constant becomes less and less important as both S^ 
and Scj, grow linearly with the scale factor. By the time 
reheating occurs, ~ ~ (2/3)A^$oa. The same rela¬ 

tion applies to perturbations in the density of nonthermal 
dark matter ^ |. 

When the Universe becomes radiation dominated, S^ 
begins to grow logarithmically and aS'^{a) is constant. If 
the transition from linear to logarithmic growth occurs 
when a = oiog, then 


S^{k > ^ rh ) 


-aiogP$o 


1 -I- In 



( 10 ) 


during radiation domination. Figures |3] and 0] show 
that this model fits the numerical solution for d;^ when 
oiog = l-29aRH, which is the same post-reheating be¬ 
havior found for nonthermal dark matter density per¬ 
turbations in Ref. [sij. When aiog is defined in this 
way, the numerical solution to Eq. m reveals that 
[kRH,h/H{a = 1)]^ = l/oiog, so Eq. (ITOl) may be rewrit¬ 
ten as 


Sx{k > fcRn) = - 



1 -kin 


/ a/ Qrh \ 

V 1-29 / ■ 
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Figure [5] shows d^ evaluated when a = lOOOaRH as a 
function of the perturbation’s wave number; we see that 
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FIG. 5: The dark matter density perturbation <5^ evaluated 
long after reheating as a function of the perturbation’s wave 
number (a = lOOOauH); modes with fc/fcuH.h < 0.001 have 
not yet entered the horizon. The solid curves show the nu¬ 
merical solutions for both values of {ov) shown in Figure [T] 
there is no distinguishable difference between these two {av) 
values. The thin dotted line shows Eq. m, while the thin 
dashed line shows Eq. (1121) with A = 9.11 and B = 0.594. 
The thick dashed line shows Eq. with the fitting functions 
A(fc/fcRH,h) and i?(fc/fcRH.h) derived for nonthermal dark mat¬ 
ter perturbations in Ref. [^ . 


Eq. (fTT|l successfully predicts the post-reheating ampli¬ 
tude of 6^ on all scales that enter the horizon during the 
EMDE (i.e. modes with k ^ A:RH,h)- Figure[5]also shows 
that the modes that have not yet entered the horizon have 
Sy, = 5/3$o, which corresponds to an adiabatic perturba¬ 
tion during radiation domination. The subhorizon modes 
in Figure [5] that entered the horizon after reheating fol¬ 
low the standard evolution of density perturbations in a 
radiation-dominated Universe with Spr ^ 5px'- 


5^{k < / crh ) = ^^>1 


A\n 



( 12 ) 


with A = 9.11, B = 0.594, and k = ahorH{ahor) [^ . 

Figure [5] indicates that the correspondence between 
perturbations in thermal and nonthermal dark matter 
extends to perturbations of all wave numbers, regardless 
of the value of (crv). Solving the perturbation equations 
for other values of (av) confirms that as long as freeze-out 
or freeze-in occurs well before reheating {Tf ^ 5Trh), the 
post-reheating evolution of matches Eq. (fTTll . There¬ 
fore, the transfer function for density perturbations de¬ 
rived in Ref. m can be applied to perturbations in ther¬ 
mal dark matter. This transfer function was derived by 
finding functions A{k/kYm,h) and B{k/kBM,h) that match 
the behavior of (^^(^/^RH.h) when inserted in Eq. (TT^ . 
The subsequent evolution of the density perturbation was 
determined by matching Eq. (1121) to the growing and de¬ 
caying modes of the Meszaros equation . It fol¬ 

lows that the impact of an EMDE on the matter transfer 


function is to take T{k) —>■ R{k)T{k), where 


R{k) 


^(fc/fcRH,h) In 


(_^\T[ R(fc/fcRH,h)aeq 
V e 3 / ahor(fe) 


9.11 In 


(^) 


0.594#^ 


(13) 


for k > 0.05fcRH,h and i? = 1 for fc < O.OS/cRH.h- In this 
expression, the subscript “eq” refers to matter-radiation 
equality, fi and /2 are constants determined by the 
baryon fraction, and a fitting function is used to find 
ahor(fc) for all modes: see Ref. [5l| for details. 

Since the reheat temperature only enters the transfer 
function through the value of fcRH.h, it is useful to derive 
an expression that directly relates fcRH,h = aRH77(aRH) 
to Trh- Numerically solving Eq. (HD for a variety 
of reheat temperatures reveals that there is a con¬ 
sistent relationship between H^ubm) and the value of 
H in a radiation-dominated universe with temperature 
Trh; 77(aRH) = 0.88y^87rG'pr-(T'RH)/3. Furthermore, 
Pr{a)a^ = 1.087pj.(TRH)aRH long after reheating, which 
quantifies the amount of entropy injected into the radia¬ 
tion bath after a = orr. It follows from the conservation 
of entropy that g^s{T)T^o? = 1.0873/'‘g*5(TRH)lRH«RH 
for all T ^ Trh, where is the effective number of de¬ 
grees of freedom that contribute to the entropy density 
s: 3*5 = s/[(27r^/45)T^]. These two results indicate that 
fcRH,h = 0.86fcRH, where /crh is the wave number of the 
mode that enters the horizon when the temperature of 
a radiation-dominated universe equals Trh. If Tq is the 
temperature of the radiation today, then 


fcRH = 


5*s(To) 


.5*5 (Trh )J Trh 


1/3 


To /87rG 


^5*RhTrjj 


(14) 



To compute the transfer function associated with a spe¬ 
cific reheat temperature, fcRH,h = 0.86A:rh(Trh) should 
be used in Eq. (USD." 

The scales that are affected by an EMDE are smaller 
than the baryon Jeans length. Since baryons do not par¬ 
ticipate in gravitational collapse on these scales, these 
perturbations grow at a slower rate than large-scale 
perturbations [^. Ref. [5l| derived a scale-dependent 
growth function D{k, z) that accounts for this suppres¬ 
sion of the growth of small-scale perturbations. This 
growth function, combined with the modified matter 
transfer function T(fc), is used to calculate the rms den- 


^ In Ref. [5ill . fcRH.h = ^RH was used to define Trh, which does 
not give the same reheat temperature as Eq. (0. The reheat 
temperatures quoted in Ref. should be divided by a factor 
of 0.86 to obtain Trh defined by Eq. 0. 
























FIG. 6: The impact of an EMDE on a{M): the present-day 
rms density perturbation in spheres containing mass M. The 
arrows show Mrh , as dehned by Eq. 113 , for the three values 
of Trh shown. In each case, the EMDE does not affect a(M) 
for M > Mrh, but it significantly enhances a{M) for smaller 
masses. 


Ref. [5l| used the Press-Schechter mass function to 
predict how this enhancement of a{M) at small masses 
would affect the microhalo population. They found that 
an EMDE dramatically enhances the number of micro¬ 
halos with masses less than Mrh- Moreover, these mi¬ 
crohalos form very early; at a redshift z = 100, over half 
of the dark matter is bound into microhalos with masses 
greater than 0.001 Mrh- The fraction of dark matter 
bound into halos is equally large at higher redshifts, but 
the microhalos are smaller. The size of the first and 
smallest microhalos is determined by the matter power 
spectrum, which usually has an exponential cut-off on 
small scales: 6^{k) oc exp[—fc^/(2A:p„t)]. For nonthermal 
dark matter, the cut-off scale kcut is determined by the 
average velocity of a dark matter particle at reheating 
and is therefore related to the mass difference between 
the scalar and dark matter [5l|. For thermal dark mat¬ 
ter, elastic scatterings between dark matter particles and 
relativistic leptons keep the dark matter and radiation in 
kinetic equilibrium after the dark matter thermally de¬ 
couples. As shown in the next section, kcut is determined 
by these interactions. 


sity perturbation in a sphere that contains a mass M: 

/ ASJc 

[D{k, z)Tik)f Pp{k)F^{kR). (15) 

In this expression, Pp{k) is the power spectrum of su¬ 
perhorizon density perturbations during radiation dom¬ 
ination, and i? = [3M/(47rpi„_o)]^^^, where pm.o is the 
present-day matter density. The window function F(kR) 
is a convolution of a top-hat window function and a Gaus¬ 
sian window function, 

3 [sin{kR) — (kR) cos{kR)] 
(fcRp ’ 

( 16 ) 

with a = 0.0001. As discussed in Ref. [^, the inclu¬ 
sion of this Gaussian window function prevents modes 
with k~^ <C R from contributing to a^{M) while leav¬ 
ing a^{M) unaltered for M > Mrh. To facilitate 
comparison with previous work, a{M) is computed us¬ 
ing the same WMAP7 [g^ cosmological parameters 
used in Ref. [5l|: — 0.135, h = 0.704 and 

A2j(fc) = (2.44 X 10-®)[fc/(0.002Mpc"^)]° ®®3-i. The 
slight difference between these parameters and the pa¬ 
rameters determined by the Planck mission [5^ would 
have a negligible impact on the results of this analysis 
and do not justify altering the transfer function. 

Figure m shows the present-day value of ct(M) for dif¬ 
ferent reheat temperatures. An EMDE alters a{M) for 
masses that are smaller than 

47r q 

/ lOMeV y/ g.sM V 10.75 

V Trh y V 10-75 ■ 


FikR) = exp 


k^{aRY 


III. KINETIC DECOUPLING AND 
FREE-STREAMING 


The momentum transfer rate for dark matter particles 
elastically scattering off relativistic particles with tem¬ 
perature T and number density rirei is 

T 

hel — (u'el'c)7Zrel 5 (1^) 

the Tjm^ factor accounts for the fact that it takes vny^jT 
elastic scatterings to significantly alter the dark matter 
particle’s momentum. As long as Pei > M, the dark mat¬ 
ter will remain in kinetic equilibrium with the radiation 
bath, and the perturbations in the dark matter density 
will be coupled to the perturbations in the radiation den¬ 
sity. Numerous studies have considered how dark mat¬ 
ter decouples from radiation during a period of radiation 
domination (g^ - Ig^ : since the radiation perturbations ex¬ 
perience damped oscillations in a radiation-dominated 
universe, dark matter perturbations on scales that en¬ 
ter the horizon prior to decoupling are suppressed. Af¬ 
ter decoupling, the residual velocities of the dark matter 
particles lead to a free-streaming effect that further sup¬ 
presses the amplitude of perturbations on scales smaller 
than the free-streaming distance. Both suppression scales 
are directly related to the temperature at which H = Pei. 
In a radiation-dominated universe, this decoupling tem¬ 
perature is defined by the relation 


rei(7kds) = 


/SttG 


/ 3 

^g»(7kds)7i;ds 


(19) 


= 32.7M 


The subscript S denotes that Ikds is the temperature 
at which the dark matter particle would decouple in a 
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standard thermal history; as shown below, the actual de¬ 
coupling temperature greatly exceeds Ik^s if dark 
matter decouples during an EMDE. However, since TkdS 
has been calculated for many dark matter candidates, 
it provides a convenient way to parametrize the micro¬ 
physics that determines the elastic scattering rate. 

If 2\dS < Tuh, then the dark matter remains cou¬ 
pled to the radiation throughout reheating. As previ¬ 
ously mentioned, the perturbations in the radiation de¬ 
crease rapidly when the Universe becomes radiation dom¬ 
inated, and if the dark matter is coupled to radiation at 
that time, its perturbations would most likely be simi¬ 
larly suppressed. However, if TkdS > 2 rh, then the dark 
matter decouples during the EMDE. When the radiation 
is at a certain temperature during the EMDE, the added 
density of the scalar field makes the expansion rate faster 
than it would be in a radiation-dominated universe with 
the same temperature. Therefore, I?(T'kds) > rei(Tkds) 
if 2kdS > Trh, and the dark matter must decouple at a 
higher temperature. 

The dependence of the kinetic decoupling temperature 
on the reheat temperature was first noted in Ref. ; if 
CTei oc as is the case for most neutralinos, they 

concluded that dark matter decouples when the radiation 
temperature is Tkd — ^kds/'^RH if TkdS > ^rh- This 
decoupling temperature follows from equating Tei to the 
expansion rate during the EMDE. The definition of orh 
implies that 


For a ^ aj, Eq. ([S]) provides an expression for the tem¬ 
perature of the radiation bath during the EMDE as a 
function of the scale factor. 


Trh 


2 g .. 


*RH 


-I 1/4 




(-) 

\aRH/ 


-3/8 


( 21 ) 


which can be used to obtain H{T) during the EMDE: 


H{T) 


5 g.{T)T* 
2 5*rh'?rh 


18nG f TT^ 

~ I ^ 


( ■57:5*rhT’j'^ 


RH 


( 22 ) 


Although H (X T'^ during the EMDE, this equation shows 
that H{T) does not equal (T/Tr,h)"^ V^87rGp7(7AH)73; 
this commonly used relation is not correct because the 
radiation temperature when a = orh does not equal 
Trh, as can be seen in Eq. (I?T|) . If (Tel oc then 

Tel oc T®, and Eq. (fT^ implies that rei(Tkd) = H{Tkd) 
during the EMDE if 


Tkd 


' g^Tkd)^ 1^^" /5 Zms_ 
.ff*(Tkds)ff*RH. V 2 Trh 


(23) 


Since Tkd > Tkds > Trh and (/*(T) monotonically de¬ 
creases as the Universe cools, we see that including 
changes in the number of relativistic degrees of freedom 



1 2 3 4 5 10 

UdsTkH 


FIG. 7: The kinetic decoupling and free-streaming scales as 
a function of Tkds/TRH, where Tkds is the kinetic decou¬ 
pling temperature in a radiation-dominated universe. The 
solid curves show fekd/^RH, where fckd is the wave number of 
the mode that enters the horizon when the expansion rate 
equals the elastic scattering rate. The long-dashed curves 
show fcfs/AiRH, where fcfg is the inverse of the free-streaming 
horizon, for the same reheat temperatures as the solid lines 
of the same color and = ISOTrh. The short-dashed and 
dotted lines show how fcfs/fcRH changes if the dark matter 
particle mass changes while Trh is held fixed at 10 GeV. 


and the exact function of H{T) during an EMDE im¬ 
plies that the dark matter decouples at an even higher 
temperature than predicted in Ref. [b^ . 

Elastic scatterings between dark matter and relativis¬ 
tic particles affect perturbation modes that enter the 
horizon prior to decoupling, so it is useful to define k^d to 
be the wave number of the mode that enters the horizon 
when T = Tkd: 

fckd _ Q(Tkd)IiI(Tkd) 

^RH Orh V^87rG'pr(TRH)/3 

5 g.jTu) f Tus V'' 

UrhI ^ * 

Figure [3 shows fckd/^RH as a function of Tkds/Tun for 
Trh = 1 GeV, 10 GeV, and 100 GeV. This ratio is 
rather insensitive to Trh at these temperatures; since 
g*{T) varies by less than 30% above temperatures of 
1 GeV, fekd/^RH — (5/2)(Tkds/TRH)®^^ when Trh > 1 
GeV. At lower reheat temperatures, however, the g* fac¬ 
tors in Eq. can increase fckd/^RH by up to a factor of 
10, so the fckd/fcRH curves shown in Figure [7] are conser¬ 
vative. Thus we see that fckd ^rh even if the standard 
decoupling temperature TkdS is only slightly higher than 
the reheat temperature. 

After the dark matter decouples from the radiation 
bath, the velocities of the dark matter particles decrease 
as 1/a. Since T oc a“®/® during the EMDE, the velocity 
of the dark matter decreases by a greater amount as the 
Universe cools if it decouples during the EMDE than if 
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it decouples in a radiation-dominated universe. The fact 
that the dark matter decouples earlier during an EMDE 
than it would in a radiation-dominated universe further 
reduces the final velocities of the dark matter particles. 
When the Universe is radiation dominated and T = Trh, 
the velocity of the dark matter particles is reduced by a 
factor of (TRH/Tkds)^^^® from what it would have been 
if the dark matter particles decoupled when T = TkdS 
in a radiation-dominated universe. Since the comoving 
free-streaming horizon is directl y re lated to the velocity 
of the dark matter particles [^1^ [t^ , 


Afs — 




da 

a^H{a) ’ 


(25) 



dark matter particles that decouple prior to reheating 
have much smaller free-streaming horizons than they 
would if they had decoupled in a radiation-dominated 
universe (b^ . 

The free-streaming horizon is evaluated by breaking 
the integral in Eq. (1^ into three parts, treating the 
transitions between the EMDE, radiation domination, 
and matter domination as instantaneous. For a < orh! 
Eq. (pni) provides H{a). After reheating, pr{a)a^ = 
1.087pi.(TRH)a^H discussion above Eq. (fHl) ]. so 

H{a) = 1.04(a/aRH)~^'\/87rGpr(TkH)/3 until matter- 
radiation equality, after which H{a) oc With these 

approximations, 


— 



Qkd 

orh/ 


1 

^RH 


X 



ln(aeq/aRH) 

1.04 



1.04 


where Okd = a(Tkd) and Ueq = 4.15 x 10“^/(Dmh^) 
is the scale factor at matter-radiation equality. Free- 
streaming suppresses perturbations with wave numbers 
greater than kfs = Using Eq. (1^ and Eq. (l23l) to 
evaluate ORn/akd and Tkd yields 


^fs 

kRH 


TOx 


pkds\ 


29/6 


U3/6 


TkdS \TrH/ 


5 5* (Tkd) 


.2 g*RH 


2/3 


L45 g.(Tkd) 


-I 1/3 


V Trh / 

^ f g^RRTRH 


1.04 


1/3 rji 


2 g*RH 

+ 1.E 


- 1 


(27) 


where r = Tkd/(T)^<js/TRH), which can be obtained from 
Eq. (l23)) . Figure [7] shows how fcfs/fcRH depends on 
itIxITbm and Tlds/TiiH- The rapid cooling of the dark 
matter particles during the EMDE makes kfs/k^H 1 
even if TkdS is just slightly larger than Trr. 

To account for the effects of free-streaming and kinetic 
decoupling, it is customary to exponentially suppress per¬ 
turbations with k > fccut, where fccut is the smaller of fckd 


FIG. 8: The rms density perturbation a{M) with 

an exponential cut-off in the matter power spectrum: 
Sm{k) oc exp[—fe^/(2fc^ut)]. The top set of curves have no cut¬ 
off, whereas the next three sets of curves have fecut/fenH = 40, 
20, and 10 as marked. The bottom set of curves shows a{M) 
in the absence of an EMDE with no cut-off; in this case Trh 
only determines Mbs- 


and fcfs While free-streaming will always sup¬ 

press perturbations, it is less clear that momentum ex¬ 
change between the dark matter particles and relativistic 
particles will have the same effect during an EMDE. As 
discussed in Section Pl B[ radiation perturbations on sub¬ 
horizon scales grow prior to reheating, so it is reasonable 
to expect that perturbations in the dark matter will ex¬ 
perience similar growth prior to decoupling. As long as 
the dark matter decouples prior to reheating, the growth 
of the matter perturbation during the EMDE will lead 
to an enhancement of the small-scale matter power spec¬ 
trum. Therefore, although modes with k > fckd > fcRH 
may not follow Eq. (ED, they will not be exponentially 
suppressed, and a more extensive analysis of kinetic de¬ 
coupling prior to reheating is required to determine their 
properties. However, Figure [7] shows that fcfs ~ fckd for 
rrix ^ 500rRH, so modes with fc > fckd are suppressed 
by free-streaming in most cases of interest. Therefore, I 
will include the effects of both free-streaming and kinetic 
decoupling by multiplying the matter transfer function 
by exp[-fcV(2fc2^t)], with the understanding that setting 
fccut equal to fckd when fckd < fcfs may underestimate the 
matter density perturbations. 

Figure [5] shows how a small-scale cut-off in the mat¬ 
ter power spectrum affects the rms density perturbation 
a { M ) for multiple values of Trh and fccut/fcRH- Insert¬ 
ing a small-scale cut-off forces a { M ) to level off for small 
masses, but it does not erase the effect of the EMDE if 
fccut/fcRH ^10. Figure [5] also shows that ^(M/Mrh) is 
not very sensitive to Trr. As seen in Figure [SI cr(M) is 
a very shallow function of M for M > Mrh, so changing 
Mrh has only a small effect on ^(M/Mrh). In Section 
IIV A[ I will show that the annihilation boost factor gen¬ 
erated by an EMDE depends only on (t(M/Mrh) and is 
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FIG. 9: The ratio of the freeze-out temperature T/ to the 
reheat temperature Trh for freeze-out scenarios that gen¬ 
erate the observed dark matter abundance. For each value 
of mx/Tnii and Trh, {(Jv) was chosen so that Eq. ([7| gives 
— 0.12. The horizontal line shows Tf/T^u = 2vT0, 
which is the minimum value of Tf/Tna required for the par¬ 
ticles to kinetically decouple after freeze-out with fecut/^RH ^ 
10. The Iow-Trh end point of each curve is set by the mini¬ 
mum value of Trh required to generate = 0.12, while the 
high- Trh end point comes from demanding that (av) < m~^. 


therefore rather insensitive to Trh- The ratio fccut/^RH 
will play a much more significant role in the calculation 
of the boost factor, however. Since both kfs/k^n and 
^kd/^RH are determined by Tkds/TRH, it is this ratio that 
will have the largest effect on the annihilation rate. 

Figure [7] shows that kcut/knu ^ 10 if Tkds/TRH ^ 2, 
which establishes a minimum value of Tkds/T rh required 
for an EMDE to enhance small-scale perturbations. If 
Tkds/TRH > 2, then Eq. (gSl) implies that Tkd/TRH ^ 
2\/T0. Since dark matter kinetically decouples after it 
thermally decouples in most cases, freeze-out scenarios 
are expected to generate enhanced small-scale perturba¬ 
tions only if T//Trh > 2-\/To. Eigure [9] shows T//Trh 
for freeze-out scenarios that generate a relic abundance of 
dark matter that matches observations and still obey the 
unitarity condition (av) < m~‘^. Since dark matter parti¬ 
cles freeze out when ^ T/10, T^/Trh ^ O.Ito^^/Trh 
in most cases. If TOj(./Trh is held fixed, there is still a mild 
dependence on Trh because as Trh increases, a slightly 
larger value of {uv) is required to generate the observed 
dark matter abundance (see Figure El), which leads to 
a slightly lower Tf/m^ ratio. Figure [H] shows that for 
= IOOTrh, it is only possible to have Tkds/TRH > 2 
and Tkd < Tf ii the reheat temperature is less than 10 
GeV, but Tkds/TRH > 2 and Tkd < T/ is possible for all 
freeze-out scenarios with mj^/TRH >150 that generate 
the observed dark matter abundance. 

Of course, the mere existence of TkdS values that sat¬ 
isfy both Tkds/TRH > 2 and Tkd < Tf does not imply 
that there are well-motivated particles that have the cross 
sections for annihilation and elastic scattering required 



FIG. 10: The evolution of the differential bound fraction given 
by Eq. (I28II for Trh = 5 GeV. At each redshift, the solid 
curve has no small-scale cut-off in the matter power spectrum, 
while the long-dashed, short-dashed, and dotted curves have 
fccut/fcRH = 40, 20, and 10, respectively. The solid vertical 
lines mark M = Mrh, and the other vertical lines show the 
cut-off masses Mcut = (47r/3)pm,ofc<lut for the curve with the 
same line type. 


to realize this scenario. Indeed, most scans of WIMP 
parameter space predict TkdS values below a few GeV 
[6^ . but these scans are limited to models that generate 
the observed abundance of dark matter in a radiation- 
dominated universe. A much smaller value of the annihi¬ 
lation cross section is required to thermally produce the 
observed abundance during an EMDE, and such parti¬ 
cles would likely have weaker elastic interactions as well. 
A search of SUSY parameter space for possible parti¬ 
cles is left for future work. In this analysis, I will re¬ 
strict myself to scenarios with T//Trh > 2-\/To and use 
the ratio /Ccut/^RH to parametrize the impact of kinetic 
decoupling and free-streaming on the dark matter per¬ 
turbations, with Eigure [7] serving as a guide for how to 
convert this ratio to a value of TkdS- 

In the absence of a small-scale cut-off in the matter 
power spectrum, small scales will be increa sing ly homo¬ 
geneous. The Press-Schechter formalism predicts 
that microhalos form when a{M) exceeds the critical lin¬ 
ear over-density 5c {5c = 1.686 for z ^ 2). Consequently, 
if a{M) increases without bound as M decreases, micro¬ 
halos will form at arbitrarily high redshift. By limiting 
the amplitude of a{M) on small scales, the A:cut/fcRH ra¬ 
tio not only determines the mass of the smallest microha- 
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M/Mrh 


demonstrated in Figure [TT] for z = 50. This insensitivity 
to Trh is not surprising given that the bound fraction is 
determined by a{M) and, as seen in Figure|Sl a{M/M yoi) 
has a very mild dependence on Trh- Therefore, while the 
reheat temperature determines the masses of the micro¬ 
halos generated by an EMDE, it does not significantly 
change the fraction of the dark matter that is contained in 
these microhalos, nor does it affect the redshift at which 
these microhalos form. The annihilation rates estimated 
in the next section depend only on the bound fraction 
and the redshift of microhalo formation, so they are not 
highly sensitive to the reheat temperature. 

IV. ANNIHILATION SIGNATURES 


FIG. 11: The differential bound fraction given by Eq. (|28l) at a 
redshift z = 50 for different values of the reheat temperature. 
As in Figure [8] the four sets of peaks correspond to different 
values of ftcut/fcan: in order of descending maxima, the curves 
show no cut-off and fecut/feaH = 40, 20, and 10. 


los, it also determines when these microhalos will form. 
Figure ITOl illustrates how /ccut/fcRH affects the evolution 
of microhalo population, which is characterized by the 
Press-Schechter differential bound fraction: 


df 

[2 

diner 


r -52 1 

d\nM 

V TT 

dlnM 

a{M, z) ^ 

2ct2(M,z)_ 


(28) 

The bound fraction = 

gives the fraction of mass that is bound in halos with 
masses between Mi and M 2 . Figure fTOl shows that most 
of the matter is already bound into microhalos at a red¬ 
shift of 400 if fccut/^RH ^ 40. For smaller values of 
kcut/knn, however, microhalo formation is postponed to 
later times: z ~ 200 for fccut = 20 A:rh and z ~ 50 
for fccut = lOfcRH- Decreasing the fccut/^RH ratio only 
slightly decreases the mass range of microhalos that are 
present at a given redshift, but if the masses of these 
microhalos are smaller than Mcut = (47i'/3)/Om^o^c)ut> 
microhalo abundance is heavily suppressed. 

Nevertheless, Figure (TU] shows that an EMDE dramati¬ 
cally enhances the microhalo abundance, even if fccut is as 
small as IO/crh- In the absence of an EMDE, the differ¬ 
ential bound fraction inherits the weak mass dependence 
of the a{M) function. Therefore, the differential bound 
fraction for M > Mrh in FigureHUl which only begins to 
noticeably deviate from zero for z ^ 25, gives an accurate 
estimate of the differential bound fraction for M ^ Mrh 
in the absence of an EMDE. The onset of standard struc¬ 
ture formation at the lower redshifts shown in Figure [TUI 
decreases the fraction of dark matter in microhalos with 
masses less than Mrh as these microhalos are absorbed 
into larger halos. 

Figure [TUI shows the bound fraction for Trh = 5 GeV, 
but the evolution of the bound fraction as a function of 
M/Mrh is the same for other reheat temperatures, as 


The annihilation rate for dark matter is 


2m2 


dV = 


{<^v) j 
2ml ■ 


(29) 


In this expression, J has been defined to characterize the 
dependence of the annihilation rate on the distribution 
of the dark matter as opposed to its microphysics. Since 
the annihilation rate is proportional to the presence 
of substructure increases the dark matter annihilation 
rate, which is often referred to as a substructure boost. 
For a halo with mass M, the boost factor 


1 + B{M) = 


J 

f p^(r) Mr^dr 


(30) 


quantifies how the annihilation rate differs from what it 
would be if the halo had a smooth density profile Pxix)- 
In the absence of an EMDE, the boost factor is often 
calculated by extrapolating the subhalo mass function 
determined by N-body simulations down to the scale 
of the smallest microhalos, as determined by the cut¬ 
off in the matter power spectrum [e.g. pTlI - lTq . Most of 
these analyses predict that B{M) is between 2 and 10 for 
M = and between 4 and 40 for M = IO^^Mq; 

Refs. obtain larger values of B because they ex¬ 

trapolate a power-law relationship between concentration 
and mass. The range in possible boost factors stems from 
uncertainties in both the subhalo mass function and the 
minimum halo mass. 

Extrapolating the subhalo mass function is only justifi¬ 
able if the matter power spectrum is featureless on scales 
between the resolution limit of the simulations and the 
cut-off scale. Since an EMDE significantly enhances the 
abundance of microhalos with Mcut < M < Mrh, it can 
generate a boost factor that greatly exceeds the value pre¬ 
dicted by extrapolating the subhalo mass function seen 
in N-body simulations. However, as seen in Figure [UJ an 
EMDE also decreases the value of (av) that generates 
the observed dark matter density. Furthermore, Figure 
[U] indicates that nix ^ IOOTrh is required to protect the 
enhanced microhalo abundance generated by an EMDE 
from the effects of dark matter free-streaming. Together, 
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Trh [GeV] 


FIG. 12: The microphysics factor that determines the an¬ 
nihilation rate ((crn)/m^) for freeze-out scenarios that gen¬ 
erate the observed dark matter abnndance. For each value 
of rUx/Tna and Trh, {(Jv) was chosen so that Eq. ([Tl gives 
= 0.12. As in Figure |9] only scenarios that obey the 
unitary condition (av) < are shown. The dotted hori¬ 

zontal lines show {av)/m^ for {av) = 3x 10“^® cm® s“®, which 
generates the observed dark matter density in the absence of 
an EMDE. 


these results indicate that the same EMDE scenarios that 
increase the boost factor will suppress (crn)/m^. In the 
absence of an EMDE, {av) = 3 x 10“^® cm® s“® generates 
the observed dark matter abundance, which implies that 

{av) _ 2.6 X 10-1® /lOOGeV 
GeV'‘ V rn^ 

Figure [T^ shows that {av)lrn^ < 10“®® for all EMDE 
scenarios with > 100. Therefore, the overall im¬ 

pact of an EMDE on the dark matter annihilation rate 
depends on the boost factor from microhalos: is it large 
enough to compensate for the reduction in {av)l'nv^ re¬ 
quired to generate the observed dark matter density? 



A. The EMDE Boost Factor 

To estimate the boost factor generated by an EMDE, 
I will assume that the microhalos that form at high red- 
shift continue to track the dark matter density as larger 
halos form. In this case, the number density of microha¬ 
los within a halo is determined by the differential bound 
fraction at high redshift: 

dUjjiicro , xPhaloW df 
d\nMf Mf d\nMf 

where Mf is the mass of the microhalo at high redshift 
{z = Zf), phaio('c) is the spherically averaged present-day 
density of the halo as a function of distance from its cen¬ 
ter, and s(r) is the fraction of microhalos that survive at 
that position within the host halo. This microhalo mass 


(32) 


function does not assume that the microhalos retain all 
of their original mass; due to tidal stripping, the present- 
day mass of a microhalo is most likely much smaller than 
Mf. Since pm,o was used to relate M and R when com¬ 
puting a{M), both Mf and phaio should include baryons. 
Although the baryons and dark matter are distributed 
very differently within the halo, this separation will not 
affect the boost factor because phaio will be integrated 
over the volume of the halo in Eq. (1351) . 

The redshift Zf should be chosen to be large enough 
that the integrated bound fraction for M < Mrh has 
not yet begun to decrease, which indicates that the mi¬ 
crohalos are being absorbed into larger halos. Figure fTOl 
indicates that Zf ^ 25 is required to characterize the 
EMDE-enhanced microhalo population prior to the on¬ 
set of standard structure formation. Beyond this lower 
limit, the best choice for Zf depends on how the evolu¬ 
tion of the bound fraction is interpreted. For example, if 
^cut/fcRH = 20, then Figure [TUI indicates that microhalos 
with M ^ O.OOIMrh begin to form at z ~ 400, and at 
z = 50, most of the dark matter is bound into microhalos 
with M ^ O.IMrh- Unfortunately, the Press-Schechter 
formalism does not tell us about the internal structure 
of the microhalos: are the microhalos present at z = 50 
smooth halos, in which case z/ ~ 50 would best char¬ 
acterize the microhalo population, or are they bound 
collections of smaller microhalos that formed at higher 
redshift? In this case, a higher value of z/ would give a 
more accurate prediction of the boost factor. Given this 
ambiguity, I will keep z/ as a free parameter until the 
final evaluation of the boost factor, at which point the 
implications of its selection will be more transparent. 

The halo’s J factor, as defined in Eq. (1^ . is the sum 
of the contribution from the EMDE-generated microhalos 
and the contribution from the rest of the host halo, which 
includes all the substructure that would be present in the 
absence of an EMDE. If /s(r) is the present-day fraction 
of dark matter that is not bound into EMDE-generated 
microhalos, 

J = >?micro + / fs (^)Px.halo 

Jo 

where Px>haio(fO is the dark matter density of the halo, 
including all the substructure that is not attributable to 
the EMDE, and 


^IdMrh / rR 




/ 


^In Mrh 


MMf) df 


/q Mf d In Mf 

fR 

X / s(r)/Ohaio(?') d®r, 

Jo 


dinM 


f 


=ZJ 


(34) 


(35) 


where Jc\{Mf) is the J factor for a single microhalo. In 
the second line, I have assumed that Jd does not de¬ 
pend on the microhalo’s position within the host halo. 
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Although the microhalo’s position will affect how much 
mass it loses due to tidal stripping, Jd is not greatly af¬ 
fected by these spatial variations because it is dominated 
by the contribution from the microhalo’s innermost re¬ 
gion. 

Numerical simulations of microhalo formation in the 
absence of an EMDE [tI, [tI, [t^ find that micro¬ 
halos with masses just above the free-streaming scale 
have extremely cuspy density profiles: p{r) oc with 
1.3 < 7 < 1.5 in the innermost regions. However, larger 
microhalos and microhalos that form in the absence of 
a cut-off in the matter power spectrum have shallower 
cusps with 7 ~ 1 in the inner region; the density profiles 
of these microhalos are well fit by a Navarro-Frenk-White 
(NEW) profile. Mergers between microhalos of similar 
mass are thought to lead to the softer NEW profile [t^. 
The evolution of the bound fraction shown in Figure [TUI 
indicates that microhalo formation after an EMDE is 
strongly hierarchical, so it is reasonable to expect that 
microhalos experience several major mergers shortly af¬ 
ter their formation. To be conservative, I adopt an NEW 
density profile for the dark matter within the microhalos: 

Px{f) = i—j —;—yr- (36) 

{r/rs){l + r/rsY 


The scale radius rg is parametrized by the halo’s con¬ 
centration: c = 7’2 oo/?'s, where r 2 oo is the radius within 
which the average matter density is 200 times the critical 
density (p 200 = 200pcrit)- The normalization of the NEW 
profile is determined by p 2 m and the concentration: 

_ r P200 _ ^ 

3 ln(l + c)-c/(l + c)’ ^ ’ 

where = Px!Pm is the dark matter fraction. If M 200 is 

the total mass enclosed within r 2 ooj this value of ps sets 
the dark matter mass enclosed within r 2 oo to be fxM 2 oo- 
Since the critical density decreases with time, r2oo is 
roughly proportional to (1 -I- z)~^ while the Universe is 
matter dominated. Since is expected to be constant, 
C200 oc (l-l-z)”^ [7l,|8^. Simulated microhalos have con¬ 
centrations between 2 and 3 at a redshift of ^30 [zHil. 
It is unknown how an EMDE would affect microhalo 
concentrations, but since the earliest forming microha¬ 
los in these simulations have the highest concentrations 
[75l |. the microhalos generated by an EMDE will proba¬ 
bly be even more concentrated. Without an EMDE, the 
first microhalos form at a redshift of 60; if concentration 
scales with formation redshift [ziill and is otherwise 
unaffected by the EMDE, then the concentration of a 
microhalo that forms at a redshift Zf will be 


C20o(2:) ~ 3 



30 \ 

1 + zJ ■ 


(38) 


Although the NEW density profile diverges at the cen¬ 
ter of the halo, the central density is limited by both 
the annihilation rate and the phase-space density, which 


cannot increase after the dark matter particles decou¬ 
ple [12 . If the annihilating dark matter at the cen¬ 
ter of the microhalo is not replenished by infalling par¬ 
ticles, the maximum density that can be achieved is 
Pmax = T^x!{{o'v)to), where to is the age of the microhalo 
[e.g. [s^. The radius at which the NEW profile exceeds 
this density is the annihilation core radius rc- If rc ^ fg 
then 



Tg mx 


(39) 


For dark matter particles that freeze out prior to 
reheating and generate the observed dark matter 
abundance, {av)jmx increases as decreases. 

If nixlT^H > 50, then {av)jmx < 8 x 10“^^ GeV“^. 
Therefore, a microhalo that formed at very high redshift 
will have 


Tg 




(I -f zfc^ 

108 ’ 


(40) 


where (1 -I- z)^c^ ^ 10® follows from Eq. (l38l) for micro¬ 
halos that form after z = 400. For microhalos with NEW 
profiles, the presence of a constant-density core decreases 
Jci by less than 10% if Vc/vg < 0.05, so the core generated 
by annihilating dark matter will not significantly affect 
the annihilation rate within a microhalo. 

For the small values of (av) /required to thermally 
generate dark matter during an EMDE, the upper limit 
on the phase-space density set by the phase-space den¬ 
sity when the dark matter kinetically decoupled imposes 
a stron ger constraint on the central density of the mi¬ 
crohalo I82 h 84]| . If dark matter kinetically decouples 
during an EMDE, the phase-space density Q ~ p/(u^)8/^ 
at decoupling is 


Qkd — Px.O^kd 
^ Px,o 


3/2 


\ 

Tkd; 

3/2 r 


T, 


RH 




3.91T8 


(41) 
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2 I 


(42) 




IO^^Mq 

Trh 

3 

Wx/Trh 

3/2 

TkdS 

pc® (km/s)® 

_80MeV_ 


50 


2'rRH 

(4 


where Tq = 2.7K is the present-day radiation tempera¬ 
ture and Eq. (lUUI) has been used to put a lower bound on 
Tkd/TRH- By Liouville’s theorem, the collisionless evolu¬ 
tion of the dark matter particles after decoupling cannot 
increase the phase-space density, so the phase-space den¬ 
sity inside the microhalo cannot exceed Qkd • The circular 
velocity Uc(r) = GM{r)/r can be used to estimate (u^) 
within the halo; if fxP — Pgi^g/r), as is the case for an 
NEW profile at radii much smaller than Vg, 


Q 


. 3/2 

Jx 


{2TrG)^/^y^Pgr^ 



(44) 
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Setting Q = Qkd defines the radius of the constant- 
density core required by Liouville’s theorem: 


— ^ 0.04 

AI/Mrh 

-IjS 

'1 -bz/' 

3/5 

TkdS 

rs 

[ 10-4 J 


[ 400 J 


.2rRH. 


- 26/5 


(45) 

for to^/Trh > 50. Figure [10] shows that most micro¬ 
halos that form after z = 400 have masses greater than 
0.00005Mr,h- Furthermore, Tkds/TRH > 2 is required 
to make A:cut/fcRH > 10, so ^ 0.05rs for all cases of 
interest. Therefore, the constant-density core required 
by Liouville’s theorem does not reduce J^\ by more than 
10%, provided that we restrict ourselves to z/ < 400. The 
core’s impact on Jd decreases even further as TO;^/Trh is 
increased beyond 50. 

The annihilation rate within a microhalo is also rather 
insensitive to the removal of mass from its outer regions; 
truncating the microhalos’ NFW density profile at re¬ 
duces Jci by only 10% if c < 3. Simulations of subha¬ 
los with NFW profiles indicate that the subhalo is com¬ 
pletely destroyed if its tidal radius is less than twice its 
scale radius 8^. Therefore, tidal stripping will most 
likely destroy a microhalo before it significantly reduces 
Jci. If the gravitational potentials of the host halo and 
the microhalo are approximated as being generated by 
point masses, the microhalo’s tidal radius is determined 
by the Roche limit: if R is the microhalo’s distance from 
the center of the host, 


Tt 


M„ 




11/3 


2Mhost(R) 


R, 


(46) 


which implies that the average density of the microhalo 
within rt is twice the average density of the host within 
R. Accounting for the extended mass profile of the host 
increases the tidal radius [s^, so the Roche limit pro¬ 
vides a simple and conservative estimate of rj. To com¬ 
pute Mhost(R), I assume that the host halo has an NFW 
profile with the concentration-mass relation provided in 
Ref. [t^ , which is based on an analysis of halo concentra¬ 
tions in N-body simulations [s^. With this host profile, 
a microhalo with c > 2 and z/ > 50 has rt > 2rs if 
i? > 1 kpc and M 200 < As the mass of the host 

halo decreases, rt/rs at a given position decreases, so 
microhalos can survive deeper within smaller hosts. For 
all host masses, rt/vs > 2 if i? > a/ i?s/(100 kpc) kpc, 
where Rs is the scale radius of the host. If the host mass 
is greater than 10^Mq, less than 1% of the host halo’s 
mass lies within a radius of i?s/(100 kpc) kpc, so the 
tidal destruction of the microhalos in the innermost re¬ 
gion of the host does not significantly affect their overall 
abundance. 

The Earth-mass microhalos that form in cosmologi¬ 
cal scenarios without an EMDE are more vulnerable to 
stripping by stellar encounters than by the host halo 
[ 7 ^ [tsI . Isvl - f^ . Although these microhalos are expected 
to lose a considerable fraction of their mass due to inter¬ 
actions with stars, some high-density cores should sur¬ 


vive [^, 111 , 193 . Eor example, simulations of microhalo- 
star scatterings indicate that the star must pass within 
0.02 pc of the microhalo’s center to significantly reduce 
its density within its scale radius and that microhalos 
should be able to survive outside the host’s inner kpc 
[78l| . Furthermore, the microhalos that are generated by 
an EMDE form far earlier than these standard microha¬ 
los, so they have much higher central densities and are 
less likely to be destroyed by stars. The microhalo sur¬ 
vival rate in the Solar neighborhood increases sharply 
with the microhalo’s mean density [^ : whereas a mi¬ 
crohalo that formed at z = 50 has only a 10% chance of 
surviving, the survival rate increases to 85% for a micro¬ 
halo that formed at z = 200. Although the impact of 
stellar encounters on EMDE-generated microhalos cer¬ 
tainly warrants further investigation, it seems reasonable 
to assume that most of these microhalos survive outside 
the innermost kiloparsec of a galactic halo. 

If we conservatively truncate the microhalo density 
profile at r = rs and neglect the slight reduction in Jd 
due to a constant-density core, 


Jc\{M,z) = 


7 23 

,M 


(47) 


= f^ — p20o{z)c{z,Zfy 


7/24 


[ln(l -f c) - c/(l -I- c)]2 ’ 
(48) 


where z/ is the formation redshift for a microhalo with 
mass M = M 2 oo{z)- (Although M 200 is not precisely the 
virial mass predicted by top-hat collapse, the difference is 
not significant, and using p 200 as the virial density eases 
comparisons with other microhalo investigations.) As a 
dark matter halo accretes mass, the product psrg remains 
roughly constant, while gro ws ap proximately logarith¬ 
mically with concentration [^, . If the EMDE micro¬ 

halos evolve in the same way prior to their absorption into 
larger halos, then Jd for a given halo will increase with 
time due to the increasing truncation radius. After the 
microhalos become subhalos, simulations of larger satel¬ 
lites indicate that rs remains roughly constant while ps 
decreases as the subhalo loses mass [85|. However, since 
the microhalos are most likely accreted relatively soon 
after their formation (while their concentrations are still 
low) most of their mass will be within 2rs and may not be 
stripped. Therefore, the possible reduction in ps due to 
stripping will be neglected, as will the possible increase 
in rs prior to the microhalo’s absorption. Jd{M) will 
be evaluated at the same high redshift as the microhalo 
mass function. 

At a given redshift, the Press-Schechter mass function 
includes halos that formed at a range of higher redshifts, 
with smaller halos forming earlier on average. There¬ 
fore, Zf and by extension c should be functions of M. 
Since c oc (1 -I- z/), increasing a halo’s formation red¬ 
shift will increase Jd{M). For example, if Jci is evalu¬ 
ated at a redshift of 50 and Eq. (1551) is used to evaluate 
the concentration, then Jd increases by a factor of 18 
as Zf increases from 50 to 400. Therefore, ignoring this 
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variation between microhalos and treating all microha¬ 
los as if they formed shortly before the bound fraction 
is evaluated provides a lower bound on the boost fac¬ 
tor. With this approximation, the concentration is inde¬ 
pendent of microhalo mass, and inserting Eq. (1481) into 
Eq. yields 


(7/72)p200(z/)cV^ 

[ln(l-I-c) - c/(l-I-c)]2 Jg d\nM 


dhiM 


l-R 

X / s(r)phaio(?’)d^r. ( 49 ) 

Jo 


I will assume that all microhalos have c = 2 at a red- 
shift Zf because this is the lowest concentration seen in 
microhalo simulations [zllzl. 

If the microhalo survival rate s(r) is taken to be 0 
inside some radius i?min and 1 outside this radius, then 


'Anicro ^ 4/^P20o(-^/)/tot(-2^/).^halo(^ ^ Amin); (^0) 

where ftot (zf) is the fraction of dark matter bound into 
microhalos with M < Mrh at a redshift Zf. Eigure fTTl 
indicates that /tot will not depend strongly on Trh for 
a hxed value of fcout/^RH- Figure [T51 shows ftot{zf) for 
reheat temperatures of 0.1 GeV, 1 GeV, and 10 GeV and 
confirms this insensitivity to Trh- Figure [T3l also shows 
how ftotizf) depends on fecut/^RH: as expected from Fig¬ 
ure [ini decreasing kcot / ^rh delays the formation of mi¬ 
crohalos and decreases their abundance. Therefore, while 
Anicro has a very mild direct dependence on the reheat 
temperature, it is very sensitive to changes in the ratio 
of the decoupling temperature to the reheating temper¬ 
ature. Since p 20 oizf) oc {1 + Zf)^, increasing the redshift 
at which ftot (zf) is nonzero will dramatically increase 
Anicro- 

Figure [13] also provides guidance regarding the opti¬ 
mal value of Zf. The lower bound on Anicro given in 
Eq. (ISOl) was derived assuming that all the microhalos 
present at z/ formed near that redshift. Since earlier- 
forming microhalos have higher densities, Eq. (|5(I)) de¬ 
pends on the product P20oizf)ftot{zf) , which is not max¬ 
imized at the same redshift that maximizes the bound 
fraction. For example, if /Ccut = 20 /crh, then only about 
5% of the dark matter is bound into microhalos at a 
redshift of 400, whereas about 75% of the dark matter 
is bound into microhalos at a redshift of 50, and yet, 
Jmicio{zf = 400) ~ 30Jmicro(^/ = 50). Therefore, if the 
microhalos that are present at z = 400 survive today, 
their contribution to J is greater than the contribution 
from the later-forming microhalos. Unfortunately, the 
fate of the earliest forming microhalos is unknown; N- 
body simulations of microhalo formation after an EMDE 
are required to determine whether they survive their ab¬ 
sorption into the larger microhalos that form at slightly 
lower redshifts. If they are destroyed, then z/ should be 
chosen to maximize ftot{zf), as this redshift gives the 
most accurate characterization of the microhalo popula¬ 
tion generated by an EMDE. Conversely, if they survive 



FIG. 13: The fraction of dark matter bound into halos with 
masses less than Mrh as a function of redshift. The top set of 
curves does not have a small-scale cut-off in the matter power 
spectrum; in this case, all of the dark matter is bound into mi¬ 
crohalos at arbitrarily high redshifts. The lower three sets of 
curves have power spectra that are exponentially suppressed 
for k > kcut- In all cases, the bound fraction decreases at low 
redshifts when some microhalos are absorbed into halos with 
masses greater than Mrh. For comparison, in the absence 
of an EMDE, the 10 “^°Mq < M < Mrh bound fraction at 
z > 50 is less than 0.05 for all possible reheat temperatures. 


as subhalos within microhalos, then zj should be chosen 
to maximize Eq. (1501) . In either case, Eq. (15(11) provides 
a lower bound on J; choosing a large value for Zf ig¬ 
nores the contribution of later-forming microhalos, while 
choosing a small value for z/ ignores the fact that some 
of the microhalos at that redshift formed much earlier 
and have higher internal densities. 

The bound fraction ftot (zf) is not directly related to 
the fraction /s{r) of dark matter that is not contained 
in microhalos today, which appears in the expression 
for J given by Eq. (l33l) . First, if all microhalos are 
destroyed within A„iin of the center of the halo, then 
fs{r < Amin) will equal one. Second, the microhalos will 
lose mass as they are tidally stripped within the larger 
halos. A halo with a concentration of 1.5 that is stripped 
to Ts retains only 61% of its original mass, and the mass 
loss increases as the concentration increases. Therefore, 
/s(r) ^ 1 — 0.6/tot(M)- If /s assumed to be constant 
for r > Amin and 1 for r < Amin, 


— Anicro 4“ Aialo(^ Amin )J~fs [Aialo Aialo ^ Amin)] , 

(51) 

where Jhaio includes the boost factor Bs{Mh) from sub¬ 
halos with M > Mrh: 


coo 

Jhaio = [1 + Bs{Mh)]47r / r^p^(r)dr 

Jo 


= [l + Bs{Mh)] 


f^MhP2oocl 


9[\n{l + Ch) - Ch/{I + Ch)Y 


(52) 

(53) 


The last line assumes that the halo has an NEW pro¬ 
file with M 200 = M/i, scale radius As, and concentra- 
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FIG. 14: The boost factor and relative boost factor generated by an EMDE. The left panel shows the boost factor as defined 
by Eq. (I3UII , while the right panel shows this boost factor divided by the boost factor from standard subhalos with masses 
greater than 10 “^^Mq [^. Both panels show the same five combinations of and /tot as labeled in the right panel. For each 
value of Zf, the largest value of /tot shown corresponds to fccut/fcRH = 40, while the second-largest values of /tot corresponds to 
fccut/fcRH = 20. Finally, the smallest boost corresponds to fecut/^RH = 10; microhalos are not present at a redshift of 400 in this 
scenario. For each value of 2 / and /tot, the dotted curve assumes that the EMDE-generated microhalos are destroyed within 
the innermost kiloparsec of the host halo, while the solid (a/ = 400) and dashed {zf — 50) curves take i?min = ^/Rs/{100 kpc) 
kpc. 


tion Ch- The concentration-mass relation provided in Ref. where 
[t^ successfully matches the concentrations of simulated 
halos with masses ranging from IO^^Mq to IO^^Mq, 
and I will continue to use it to evaluate the present- 
day concentration of the host halo. Since subhalos with 
M > Mrh are less dense than the EMDE-generated mi- 
crohalos, they should not survive at radii smaller than 
Rmiii) which implies that 


J{y) = 


-h 3y 
3(?/-h 1)3 


(56) 


pRmin 

^haio(r<Rniin) = dTT / r'^p‘^{r)dr, (54) 

Jo 

_ _ f^MhP200cl _ / Rmin \ 

3[\n{l + Ch)-Ch/{l + Ch)]'^ V i?s y ’ 

(55) 


The total boost factor generated by an EMDE is ob¬ 
tained by dividing Eq. dnH) by 47r /“ r^p^(r)dr: 


l + B^36/tot(^/) 


P200{zf) 

' Mh{r> Rniin)' 

[ln(l -b Ch) - Ch/{1 + Ch)Y 

P200 

[ Mh \ 

c3 


[1 + B,{MnM + (^) (1 - fs)- 

^^ 


Figure [14] shows 1-1-R as a function of Mh, along with the 
relative boost from an EMDE: (1 -l-i?)/(l -t-i?s). In both 
cases, Bs(Mh) was evaluated using the fitting function in 
Ref. [7^ for the boost factor from subhalos with masses 
larger than 1O“®M0 and a mass function dN/dm oc m~^ 
multiplied by a factor of 2.5 to match their results for the 
boost factor with a minimum subhalo mass of 10 
As seen in Figure 1121 the most promising scenarios for 
detection have reheat temperatures between 350 MeV 
and 10 GeV, which implies that Mrh is between IO^^Mq 


and Since B{M) changes by only a factor of 

2.5 when the minimum subhalo mass increases by six 
orders of magnitude, using M^in = instead of 

A/min = A7 rh will not significantly affect the boost factor 
for these reheat temperatures. In Figure (Td] fs has been 
set to 1 —0.6/tot(-Z/), but since the Eq. ((57)) is dominated 
by the first term for the cases shown, this choice has very 
little impact. 

In Figure [Til Zf and /tot are chosen to correspond to 
optimistic and pessimistic scenarios for fccut/^RH = 40 
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and fecut/^RH = 20. (If there is no small-scale cut-off 
in the matter power spectrum, an EMDE will generate 
microhalos at arbitrarily high redshifts, which will lead 
to arbitrarily large boost factors.) In the optimistic case, 
Zf = 400 is chosen to maximize Eq. (1351) . while in the 
pessimistic case, the microhalos present at high redshift 
are assumed to be destroyed and Zf is taken to be 50. 
Figure [m also shows Zf = 50 for fccut/^RH = 10; in 
this case, choosing a higher Zf does not significantly in¬ 
crease the boost factor. Figure [T3] shows that larger halos 
have larger absolute boost factors, but the relative boost 
from an EMDE is nearly independent of mass if i?niin is 
small enough that Mh{r > i?min) — Mh- The origin of 
this behavior can be seen in Eq. m-- larger halos have 
lower concentrations than smaller halos, which means 
that high-density microhalos have a larger impact on 
the dark matter annihilation rate. However, when stan¬ 
dard subhalos are included, they dominate Jhaio- Conse¬ 
quently, the relative boost from an EMDE is attributable 
to the microhalos’ enhancement of the annihilation rate 
within the subhalos, which does not depend on the mass 
of the host halo. 


B. Detection outlook 

For dark matter particles with masses less than a 
TeV, the most stringent constraints on the dark mat¬ 
ter annihilation rate come from observations of dwarf 
spheroidal galaxies (dSphs) by the Fermi Large Area 
Telescope (Fermi-LAT). The Fermi-LAT Collaboration 
has recently updated these constraints using six years of 
observations of fifteen dSphs [s^. For dark matter an¬ 
nihilating into b quarks, (av) > 2 x 10“^®cm^s“^ is ex¬ 
cluded at 95% confidence level (CL) for = 100 GeV, 
and {av) > 2 x 10“^®cm^s“^ is excluded at 95% CL 
for = 1 TeV. Slightly weaker constraints are found 
if dark matter annihilates into r leptons. For particles 
heavier than 1 TeV, air Cherenkov telescopes provide 
more stringent constraints on (av)- For = 10 TeV, 
MAGIC observations of Segue 1 [9^ and H.E.S.S. obser¬ 
vations of five dSphs [9^ exclude (av) > 10“^^cm^s“^ 
at 95% CL if dark matter annihilates into t leptons, and 
observations of the Galactic center by H.E.S.S. exclude 
(av) > 10“^^cm^s“^ at 95% CL if dark matter anni¬ 
hilates into b quarks [o^- However, these H.E.S.S. ob¬ 
servations are conhned to a region within 150 pc of the 
Galactic center, which is expected to be devoid of mi¬ 
crohalos. Based solely on observations of dSphs, Fermi- 
LAT provides the strongest constraints on heavy dark 
matter particles annihilating into b quarks, excluding 
((to ) > 4 X 10“^"'^ cm^ s“^ at 95% GL for = 10 TeV 

None of these constraints on (av) from dSphs include 
a boost factor from substructure; given the uncertainty 
in the boost factor, it is considered more conservative to 
omit it entirely. Unfortunately, the boost factors shown 
in Figure [T3] cannot be readily applied to observations 


of dSphs because these observations are confined to the 
dSphs’ central regions. The Fermi-LAT analysis only 
considered emission from within 0.5° of the dSphs’ cen¬ 
ters, which corresponds to a radius of 200 pc for Segue 
1 and 280 pc for Ursa Major, the two nearest dSphs. If 
these galaxies have NFW profiles with the concentration- 
mass relation given by Ref. [t^, then microhalos are 
expected to survive at radii greater than 30 or 40 par¬ 
secs [i?min = \/i?s/(100kpc) kpc] in these systems, so 
an EMDE will enhance the annihilation rate within this 
region. However, Jmicro is proportional to the halo mass, 
and only 15% of the mass of Ursa Major and 25% of the 
mass of Segue 1 lies within 0.5° of their centers. Mean¬ 
while, restricting Jhaio to this region reduces Jhaio by less 
than 7%. Therefore, Fermi-LAT’s assumption of a lim¬ 
ited annihilation region will reduce the microhalo boost 
factor by the fraction of mass that is included, and the 
boost factors shown in Figure [T3] should be reduced by 
roughly a factor of 10 when applied to dSphs. 


This reduction in the boost factor for dSph obser¬ 
vations makes it worthwhile to consider the nominally 
weaker constraints from Fermi-LAT’s observations of the 
isotropic gamma-ray background (IGRB) [s^. Conserva¬ 
tive constraints are derived by assuming that the IGRB 
originates solely from dark matter annihilations. The 
95% GL upper limits on (av) for annihilation into b 
quarks are 7 x 10“^^cm^s“^, 3 x 10“^“^ cm^ s“^, and 
5 X 10“^"'’cm^ s“^ for = 100 GeV, 1 TeV and 10 TeV, 
respectively. The upper limits for annihilation to t lep¬ 
tons differ by less than an order of magnitude. To account 
for substructure. Ref. uses the same Bs{M) function 
from Ref. [t^ as I do, but with Mmin = 10~^h~^MQ. 
Since the annihilation signal is proportional to the boost 
factor, decreasing Mmin to would decrease all 

these upper limits on (av) by a factor of ^2/5. Ref. 
[ 5 ^ also provides an estimate of Fermi-LAT’s sensitivity 
reach by constraining a possible dark matter annihilation 
signal that lies on top of a simple model for the contribu¬ 
tion to the IGRB from other astrophysical sources. These 
sensitivity bounds forecast the constraints on (av) that 
could be derived from the IGRB if the astrophysical back¬ 
ground were fully understood. For ^ 1 TeV, these 
upper limits are about an order of magnitude stronger 
than the conservative upper bounds on (av). 


To use these limits on the dark matter annihilation 
rate to constrain EMDE cosmologies, we must deter¬ 
mine the appropriate boost factor. Since no boost fac¬ 
tor was included in the analysis of the dSphs, the ratio 
{av)jm‘^ should be multiplied by the total boost factor 
1 + B prior to comparison with the derived upper lim¬ 
its on {av)lm'^. The dSphs with the largest J factors 
have masses around 10^Mq. As seen in the left panel 
of Figure [HI 1 + B ^ 200,000 for a halo with a mass of 
10®M 0 . However, as discussed earlier, the Fermi-LAT’s 
limited annihilation region effectively reduces the boost 
factor by roughly a factor of 10. Therefore, the effective 
boost factor for dSphs is about 20,000 for Zf = 400 and 
ftot{zf) = 0.6, which corresponds to fecut/^RH = 40, as 
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FIG. 15: The ratio {av)lm^ times the relative boost fac¬ 
tor from EMDE-generated microhalos, along with current up¬ 
per bounds on {av)lm'^ from Fermi-LAT observations. The 
curves show values of (1 -I- B){av)/m\ and Trh that gener¬ 
ate the observed dark matter abundance for two values of 
mj(./rRH and two values of the boost factor B. The solid 
curves have 1 + B = 20,000 (appropriate for dSphs), and 
the dashed curves have 1 + B = 75,000 (appropriate for the 
IGRB). The line segments show existing upper bounds on 
{av)/m'^ for = 100 GeV, 1 TeV, and 10 TeV from Fermi- 
LAT observations of dSphs and the IGRB [s^. The dot¬ 
ted segments labeled IGRB (S) show the upper bounds on 
((Tii)/m^ that could be obtained if the astrophysical contri¬ 
butions to the IGRB can be removed [13] . These line segments 
are positioned so that = IOOTrh; they may be directly 
compared to the = IOOTrh curves. 

seen in Figure fCT 

The solid curves in Figure ITS] show 20,000 ((T'(;)/to^ for 
EMDE scenarios that generate the observed dark matter 
relic abundance along with the upper bounds on (av) / 
for = 100 GeV, 1 TeV, and 10 TeV derived from 
Fermi-LAT observations of dSphs. Figure [12] indicates 
that current limits on gamma-ray emission from dSphs 
are already capable of ruling out EMDE scenarios with 
^ IOOTrh, Trh ^ 10 GeV, and fccut/^RH ^ 40, 
which corresponds to TkdS ^ 3 Trh. Although a more 
detailed computation of the boost factor for each indi¬ 
vidual dSph would be required to firmly establish these 
constraints, this preliminary estimate encourages such an 
undertaking. 

Figure[T5|also illustrates the potential power of the the 
IGRB to constrain EMDE scenarios. Since the Fermi- 
LAT analysis of the IGRB included a boost factor from 
standard subhalos with masses greater than 
the ratio {av)(m'^ should be multiplied by the relative 
boost factor 

1 + B ^5 1 + B 

1 + BsiM^in = IO- 6 M 0 ) “21 + BsiM^in = 10 - 12 Mq) ■ 

(58) 

If Zf = 400 and ftot{zf) = 0.6, then the right panel 
Figure [TH tells us that the appropriate boost factor is 
5/2x30,000=75,000. In Figure |T21 this boost factor cor¬ 


responds to the dashed curves. Although 75,000(crr!)/TO^ 
falls below the current IGRB bounds on the annihilation 
rate for all interesting values of the IGRB may 

be able to constrain scenarios with fecut/^RH ^ 40 be¬ 
cause microhalos form at redshifts higher than 400 in 
these scenarios. The boost factor may be significantly 
greater than 75,000 if these early-forming microhalos sur¬ 
vive as subhalos within the microhalos present at a red- 
shift of 400. Even if the microhalos at z/ = 400 are 
devoid of substructure, accounting for the earlier forma¬ 
tion times of the smaller microhalos present at that red- 
shift may increase the boost factor by the factor of two 
required to reach the observational constraints. Figure 
US) also shows that the Fermi-LAT observations of the 
IGRB are easily sensitive enough to detect gamma rays 
from dark matter annihilations in EMDE scenarios with 
~ IOOTrh if the astrophysical background can be 
removed. 

Since the ratio (av) for scenarios that generate the 
observed relic abundance decreases rapidly as 
increases (see Figure [T^ . it is unlikely that gamma-ray 
observations will be able to constrain EMDE scenarios 
with m^/TRH 200. As shown in Figure [T21 even 
m^/TRH = 150 lies well within current observational 
bounds on the dark matter annihilation rate. Further¬ 
more, the large boost factors (1 + R ^ 20,000) required 
to bring the annihilation rate above current observational 
bounds cannot be obtained if z/ ^ 100. Since microha¬ 
los do not form before z ~ 100 if fccut/fcRH = 10, as 
shown in Figure [T31 only scenarios with /Ccut/^RH 20 
can generate large enough boost factors to saturate cur¬ 
rent observational constraints. 

Although the estimate of the boost factor developed in 
this section indicates that fecut/^RH 40 is required to 
obtain boost factors greater than 10^, it is important to 
remember that the boost factor may greatly exceed this 
estimate. It is possible that the boost from later-forming 
microhalos, which is neglected when we take z/ = 400, 
is sufficient to make 1 + R ~ 10^ with fecut/^RH — 20. 
In light of this uncertainty, current observations may be 
able to constrain EMDE scenarios with fecut/^RH ^ 20, 
which corresponds to TkdS 2.1Trh for ~ IOOTrh. 

Any improvement in the observational constraints on 
the dark matter annihilation rate could extend the reach 
of these observations to EMDE scenarios with larger 
W;^/Trh ratios or lower decoupling temperatures. Four 
more years of Fermi-LAT observations of the dSph sam¬ 
ple used in Ref. [s^ is expected to moderately strengthen 
the upper bounds on {av)lm‘^ [o^. Including the ultra¬ 
faint dSph candidates that were recently discovered in 
the Dark Energy Survey [H, llOdj] could provide a more 
dramatic improvement when t he d ark matter content 
of these systems is established |lQl| . Finally, a better 
understanding of the astrophysical contributions to the 
IGRB would enable their removal, which would make the 
IGRB limits on dark matter annihilation far more sensi¬ 
tive, as seen in Figure [T21 [ 5 ^. 
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V. SUMMARY AND DISCUSSION 

A transient period of effective matter domination prior 
to BBN is a generic feature of inflation theories and 
string theories that include gravitationally coupled mod¬ 
uli fields. Moreover, the lightest modulus is expected to 
dominate the energy density of the Universe at tempera¬ 
tures higher than a few hundred GeV if supersymmetry 
provides even a partial solution to the electroweak hier¬ 
archy problem [e.g. [l^. This delayed onset of radiation 
domination dramatically alters the relationship between 
the properties of the dark matter particle and its present- 
day density [29l - [3^ liol - ld^ . Consequently, limits on the 
dark matter annihilation cross section cannot rule out a 
thermal origin for dark matter as long as the reheat tem¬ 
perature is unknown. Unfortunately, the period between 
inflation and BBN is difficult to probe; nonthermal ex¬ 
pansion histories are fully consistent with the primordial 
abundance of light elements and cosmological observa¬ 
tions of large-scale structure provided that the reheat 
temperature is at least 3 MeV M- 

Microhalos offer a new window into the early Universe 
because they form from density perturbations that enter 
the cosmological horizon prior to BBN. Earlier analyses 
of perturbation evolution during reheating [H, Is^l 
have established that dark matter density perturbations 
grow linearly with the scale factor prior to reheating if 
the dark matter is a decay product of the scalar field 
that dominates the Universe during the EMDE. In this 
work, I have extended this investigation to include dark 
matter that is thermally produced during the EMDE. Al¬ 
though thermal and nonthermal dark matter perturba¬ 
tions evolve differently at early times, the two scenarios 
converge after dark matter thermally decouples. Fur¬ 
thermore, the final amplitude of the perturbations does 
not depend on whether or not the dark matter reaches 
thermal equilibrium prior to decoupling. Consequently, 
the post-reheating matter power spectrum for both the 
freeze-in and freeze-out production mechanisms matches 
the power spectrum derived for nonannihilating nonther¬ 
mal dark matter in Ref. M- In all cases, the linear 
growth of matter perturbations during the EMDE signif¬ 
icantly enhances the amplitude of small-scale dark matter 
density perturbations. 

The velocity dispersion of the dark matter particles 
and their elastic scatterings with relativistic particles can 
suppress small-scale perturbations and erase the pertur¬ 
bations that grow during an EMDE. This is a crucial con¬ 
cern for nonthermal dark matter because free-streaming 
erases the perturbations that grow during an EMDE if 
the dark matter particles are relativistic at reheating 
pm. If dark matter is produced thermally, however, 
it may kinetically decouple during the EMDE. In this 
case, the cut-off in the matter power spectrum from free- 
streaming and elastic scatterings does not suppress all 
the scales that are affected by an EMDE. Furthermore, 
dark matter particles decouple earlier and cool faster dur¬ 
ing an EMDE than they would in a radiation-dominated 


universe, so a small relative difference between the stan¬ 
dard decoupling temperature and the reheat tempera¬ 
ture leads to a much larger relative difference between 
the cut-off scale and the horizon size at reheating. If 
dark matter kinetically decouples at a temperature TkdS 
in a radiation-dominated universe, then TkdS ^ STrh is 
sufficient to preserve the enhanced inhomogeneity gener¬ 
ated during an EMDE. Consequently, an EMDE signif¬ 
icantly enhances the microhalo population predicted by 
the Press-Schechter mass function if dark matter ther¬ 
mally and kinetically decouples prior to reheating. For 
TkdS ^ 2rRH, most of the dark matter is bound into 
microhalos at a redshift of 100. These early-forming mi¬ 
crohalos are far denser than the microhalos that form in 
purely thermal cosmologies, so they should survive out¬ 
side the innermost region of a larger halo. 

A thermal origin for dark matter provides a means 
to detect these microhalos: dark matter annihilations 
within their dense cores will boost the overall annihi¬ 
lation rate within their host halos. Unfortunately, this 
substructure does not automatically lead to an overall 
enhancement of the annihilation rate because a smaller 
value of the annihilation cross section is required to gen¬ 
erate the observed dark matter abundance during an 
EMDE. If the dark matter particle freezes out during 
an EMDE and ^ IOOTrh, the ratio {av)fm'^ re¬ 
quired to generate the observed dark matter density is at 
least three orders of magnitude smaller than it is for a 
100 GeV particle that freezes out during radiation dom¬ 
ination. Freeze-in scenarios demand even smaller anni¬ 
hilation cross sections for the same particle masses and 
reheat temperatures. The resulting reduction in the dark 
matter annihilation rate is so severe that annihilations of 
thermal relics generated during an EMDE were thought 
to be undetectable, but an EMDE’s effect on the micro¬ 
halo population could bring these scenarios within reach 
of current gamma-ray observations. 

To estimate the substructure boost factor from an 
EMDE, Press-Schechter mass functions were used to pre¬ 
dict the fraction of dark matter that is contained in mi¬ 
crohalos at a certain redshift. All the microhalos at that 
redshift were assumed to have NEW profiles with c = 2. 
The density profile within the scale radius was assumed 
to be unaltered by the microhalos’ subsequent absorption 
into their host halos, except in the innermost region of the 
host halo where the microhalos’ tidal radii are less than 
twice their scale radii. The microhalos within this region 
are most likely destroyed, so they do not contribute to 
the boost factor. Outside this region, the distribution 
of the microhalos was assumed to follow the dark mat¬ 
ter density profile of the host halo, which was taken to 
be an NEW profile with the concentration-mass relation 
proposed in Ref. [t^ . 

The resulting boost factor is largely insensitive to the 
reheat temperature, but it depends very strongly on the 
redshift at which the microhalo mass function is evalu¬ 
ated. If 90% of the dark matter is contained in microhalos 
at a redshift of 50 (as is predicted for TkdS ^ STrh), the 
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relative boost factor from an EMDE is 100, which is not 
large enough to overcome the suppression of the annihi¬ 
lation cross section. The evolution of the Press-Schechter 
mass function following an EMDE is strongly hierarchi¬ 
cal, however, and microhalos are also common at very 
high redshift {z ^ 400 for TkdS ^ 3 Trh)- If these early- 
forming microhalos survive as subhalos within larger mi- 
crohalos, then the relative boost factor from an EMDE 
can exceed 30,000. This boost factor is sufficient to make 
the dark matter annihilation rate within dSphs exceed 
the limits established by Eermi-LAT if ~ IOOTrh 
and Trh ^10 GeV. It would be interesting to search for 
dark matter candidates that have the masses and cross 
sections required to realize this scenario, as these models 
may already be in tension with observations. The IGRB 
[STjl may also constrain these scenarios; increasing the 
estimated boost factor by a factor of two would be suffi¬ 
cient to saturate current IGRB bounds on the dark mat¬ 
ter annihilation rate. Improving the characterization of 
the microhalo population could easily increase the boost 
factor by this amount. 

The primary source of uncertainty in the EMDE boost 
factor is the internal structure and substructure content 
of the microhalos that form from the enhanced small- 
scale perturbations. The Press-Schechter formalism only 
provides a prediction for the mass function; it cannot de¬ 
termine if the numerous microhalos present at a redshift 
of 400 survive within the larger microhalos that contain 
most of the dark matter at a redshift of 50. Numer¬ 
ical simulations of microhalo formation in EMDE cos¬ 
mologies, similar to the simulations of microhalo forma¬ 
tion that have already been done for thermal histories 
[zl Izi, Izllzi, are required to determine the internal 
structure of the microhalos. The extremely hierarchical 
nature of microhalo formation following an EMDE facili¬ 
tates these simulations because EMDE-generated micro¬ 
halos form much earlier than significantly larger halos. 
Therefore, it should be possible to extract a microhalo 
mass function from simulations in EMDE cosmologies, 
providing a way to check the predictions of the Press- 
Schechter formalism. 

The annihilation boost factor from an EMDE is also 
strongly dependent on the dark matter decoupling scale; 
a cut-off in the small-scale power spectrum determines 
when the first microhalos form. In EMDE cosmologies 
without a small-scale cut-off, all of the dark matter is 
bound into microhalos at arbitrarily high redshift, which 
could lead to very high annihilation rates. In this work, 
the small-scale power spectrum was assumed to be ex¬ 
ponentially suppressed on scales smaller than the free- 
streaming length and the horizon size at kinetic decou¬ 
pling. This treatment is based on analyses of dark matter 
decoupling during radiation domination [e.g. How¬ 
ever, the radiation perturbation evolves very differently 
during an EMDE , so it is possible that elastic scatter¬ 
ings between dark matter particles and relativistic parti¬ 
cles do not have the same effect on perturbations as they 
do during radiation domination. An analysis of pertur¬ 


bation evolu tion through kinetic decoupling during an 
EMDE |l02j | will determine the exact relationship be¬ 
tween the cut-off scale and the elastic scattering cross 
section and may extend the reach of gamma-ray obser¬ 
vations beyond the TkdS ^ 2 Trh limit established here. 

If gamma rays from dark matter annihilations are de¬ 
tected, we must disentangle the dark matter particle’s 
properties from the potential boost factor from an en¬ 
hanced microhalo population. Eortunately, the impact of 
an EMDE on the spatial variation of the dark matter an¬ 
nihilation rate cannot be mimicked by simply increasing 
the annihilation cross section. Since microhalos are not 
expected to survive within the innermost kiloparsec of 
the Galaxy, emission from the Galactic center will not be 
affected by an EMDE. Therefore, a distinctive signature 
of an EMDE would be an annihilation signal observed 
in both dSphs and the IGRB with no corresponding sig¬ 
nal from the Galactic center. Furthermore, if the boost 
factor from an EMDE is significant, then the gamma- 
ray signal from dark matter annihilations within halos 
will be proportional to the number density of microha¬ 
los, which will scale with the density of the halo and not 
its square. As a result, an EMDE can make the gamma- 
ray emission from dark matter annihilations resemble the 
emission from decaying dark matter outside of the halo’s 
innermost region, but the strength of the signal from the 
halo center would be inconsistent with this interpreta¬ 
tion. Finally, an enhanced microhalo population would 
also affect the gamma-ray angular power spectrum [l03j| , 
which may provide an additional way to determine if the 
annihilation rate has been boosted by an EMDE. 

An EMDE widens the field of dark matter candi¬ 
dates; particles with small annihilation cross sections 
that would be overproduced in a radiation-dominated 
universe are viable dark matter candidates if they freeze 
out before reheating. I have shown that the effect of an 
EMDE on small-scale perturbations in these scenarios 
provides an additional observational signature that can 
be used to constrain these models. If dark matter ki- 
netically decouples before reheating, the growth of per¬ 
turbations during the EMDE leads to an abundance of 
microhalos that boosts the dark matter annihilation rate 
and alters its spatial variation. Therefore, both the ab¬ 
sence of gamma rays from dark matter annihilation and 
the properties of a signal detected in the future can pro¬ 
vide a new window into the early Universe, increasing 
our understanding of inflation, reheating, and the origins 
of dark matter. 
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Appendix A: Derivation of the Perturbation 
Equations 


The equations that govern the evolution of the per¬ 
turbations are derived by perturbing the covariant form 
of th e energy-transfer equations given in Eq. ([IJ. Ref. 


104j used this method to derive equations for the evolu¬ 


tion of density perturbations in annihilating dark matter 
on superhorizon scales during reheating in the curvaton 
scenario. The same approach was used in Ref. m to 
obtain the evolution equation for density and velocity 
perturbations on all scales for nonthermal dark matter. 
Refs. recently extended the analysis of Ref. [5l| 

to include dark matter annihilations. In this appendix, I 
review this derivation and apply it to the case that dark 
matter is not generated during scalar decays. 

The oscillating scalar field, the radiation, and the 
dark matter are all treated as perfect fluids with energy- 
momentum tensors 


T>^'' = {p+p)uV +pg^^'', (Al) 


where p and p are the fluid’s density and pressure, respec¬ 
tively, and = dx'^/dA is its four-velocity. The dark 
matter and the oscillating scalar fields are both treated 
as pressureless fluids, while the radiation has p = p/3. 
Since the scalar field decays into radiation and the dark 
matter can self-annihilate into relativistic particles, these 
three fluids exchange energy, as described in Eq. ([T]) . This 
energy exchange can be expressed covariantly as 

V^(Wr^)=QW, (A2) 

where i denotes the individual fluids. In the absence of 
spatial variations, 

V;. = -p^ - 5H{p, + Pi), (A3) 

= 0, (A4) 

where a dot denotes differentiation with respect to proper 
time. It is also useful to note that = —u^p and 

TvxT^^up = Uyp^. If the fields are homogeneous, then 
Mo = —1 and Ui = 0. 

It follows from Eq. o that 


= T^m^T^ = -M^pT^, (A5a) 

Qij) =-qW+L,,, (A5b) 

= -T., (A5c) 


where 

L,= 


(ctm) 


rriy 


{x)t X _(x,eq) (x,eq)/j^A/3 


/ 2 2 


rriy 


(Px - 


(x) 


,(x) 


(A6) 


captures the energy exchanged between the radiation 
bath and the dark matter through annihilation and pair 
production. In this three-fluid model of reheating, 

+ (A7) 


as required by the conservation of energy and momen¬ 
tum. 

The perturbation equations are obtained by evaluating 
Eq. (IA2I) with the perturbed FRW metric 

ds^ = —(1 -I- 24')dt^ -I- a^(t)Sij{l + 24>)da;Ma;'’ (A8) 


and perturbations in the density of each fluid: 
Pi{t,x) = Pi{t)[l + 6i{t,x)]. I also introduce perturba¬ 
tions to the four-velocity of each fluid: Mq = —(1 -I- 4') 
and Ujpi) = a^SkjV^iy where = dx^/dt is the peculiar 
fluid velocity of the fth fluid in comoving coordinates. It 
follows that 

Q[,'^)=r^pO(l + <5^ + 4/), (A9) 

Qf = -T^p%aHujvl, (AlO) 

to first order in the perturbations. Note that is 

first order in the perturbations, while ^ has both a 
zeroth-order component = T^p^] and a first- 

order component = T^p'^{S^ + 4')]. In addition, 

L, = -^Hpir(l + 2S, + <l) 

rri^ 

“ (7'x,eq)^(l + 2^X.eq + ^)] I (^H) 

'"■X 

where d^,eq is the perturbation in the equilibrium density 

of the dark matter defined in Eq. (jH). Like Q^it\ Lj is 
a first-order quantity, while Lq has both a zeroth-order 
component lI°^ = -{{av)/m^) [(p°)^ - (p°,eq)^]> and a 
first-order component 


4'^ = [(Px)'(24 + 'f') - (Px.eq)"(2<Jx,eq + 'I')] . 

(A13) 

The V = Q component of Eq. (IA2|) implies that each 
fluid obeys the equation 


d(5 9 d4> 1 

—-f(l-bM;)-+3(I-bM;)— = — 
at a at 




, (A14) 


where w = p/p is the fluid’s equation of state parameter, 
9 = a diV^ is the divergence of the fluid’s physical velocity, 
and and are the zeroth-order and first-order 
components of Qo for tiiis fluid. The divergence of the 
spatial components of Eq. (IA2I) implies that 


w 


V^S 


d9 

- + (l-3w)II9 + -+ —- 

dr a 1 -I- w a 

1 


diQi 

i(l -I- w) 




(A15) 
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where 

a(l + w) 
diLi 

a(l + w) 


{av) 6^ 


my 1 + w 


(1 + w) 




A T 

l + w^° 


(0) 


(A16) 




(A17) 


Applying Eqs. (jA14l) and (IA15I) to the scalar, radia¬ 
tion, and dark matter perturbations yields 


J 


a at 


dt 


d6>0 

dt 


H0, 




= 0, 


a 


dt 


A +3:iA =_ 

a dt 'mxPx 


[-^{(p;)" - - (pifSx + - <«)]. 


d(5r 4 6»^ d$ r. pitX X , .T,^ , W 

17 + 31- + “* + + 


(A18a) 

(A18b) 

(A18c) 

(A18d) 


dOr 

dt 


1 V^Sr 


4 a 


= (-e^-Or 1 - 


{<^v) 

mxP^r 


MiP^xf - K.eq)"} + (Px)'(2<^X - - (Px.eq)"(2<5x.eq ” <5.)] , 

(A18e) 

(A18f) 


[(Px)" - (/’X.eq)"] [<^r - \e^ . 


The Einstein equations for the gravitational potentials dt 
and $ imply 

^ + 3i7 (^Hd! - ^) = -47rG {p^S^ + p^ + p^Sy) . 

(A19) 

Since the radiation fluid is tightly coupled, the 
anisotropic stress is negligible, and the Einstein equa¬ 
tions demand that $ = — 


Appendix B: Perturbation Initial Conditions 

The solution to Eq. (IA18I) is obtained by expressing 
the variables as functions of a instead of t and then 
numerically integrating the resulting set of coupled dif¬ 
ferential equations for a single plane-wave perturbation 
mode with wave number k from a = 1 to some value of 
a > orh = where Hi = H{a is the 

initial value of the Hubble parameter (Sf], |H, [sj. The 
integration begins when the mode is outside the cosmo¬ 
logical horizon: k <C aHi. The addition of dark matter 
annihilations does not affect the evolution of the pertur¬ 
bations in the scalar field, so the evolution of the scalar 
perturbations and the gravitational potential $ during 
the EMDE is the same that as derived in Ref. [Slj: 

$ = $ 0 , (Bla) 


6^ = 24)o -I- ^P$oa, (Bib) 

^ = -‘^k'^^ox/a, (Blc) 

where k = k/Hi. Since the number of relativistic parti¬ 

cles created in scalar decays greatly exceeds the number 
of relativistic particles created or destroyed via dark mat¬ 
ter annihilations, adding dark matter annihilations also 
does not affect the evolution of the perturbations in the 
radiation. For superhorizon modes |5l| . 

(5r- = $0 + (B2a) 

63 

^ = (B2b) 

111 o 

The background equations given by Eq. © apply to 
both relativistic and nonrelativistic dark matter, but the 
perturbation equations given by Eq. (jA18l) are only ap¬ 
plicable to nonrelativistic dark matter. Therefore, the 
numerical integration of Eq. (IA18I) must start at a tem¬ 
perature such that T/my < 1. For the freeze-out sce¬ 
nario, however, we can only capture the evolution of the 
perturbations during freeze-out if the initial temperature 
exceeds the freeze-out temperature. Since Tf ^ rnyjb, 
I start the integration of Eq. (jA18l) when m^/T = 3 for 
values of (av) that allow the dark matter to reach thermal 
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equilibrium. While the dark matter is in thermal equi¬ 
librium, the individual terms on the rhs of Eq. (lAlScI) 
are much larger than the terms on the Ihs, so the rhs 
terms must nearly sum to zero. If ~ Px,ev the rhs of 
Eq. llAlScI) vanishes only if 

~ (B3a) 


which establishes the initial condition for 5^. Further¬ 
more, this suite of initial conditions for 5^, <5^, and 
forms set of adiabatic perturbations because H{pi/pi)5i 
is the same for all three components. 

It is less clear what the appropriate initial condi¬ 
tion for (5^ is if the dark matter never reaches thermal 
equilibrium. For these values of (cm), the value of 5^ 
on superhorizon scales is chosen to make the perturba¬ 
tions adiabatic. Equations (jBlIl and (IB2I) already ensure 
that HS^{pff,/p^) = H6r{pr/ Pr) during the EMDE, since 
/90 oc a~^ and pr oc prior to reheating. The dark 

matter perturbation preserves adiabaticity if 




I^ ap'^ia) 
3^ Pxia) ’ 


(B3b) 


where p'^ia) is obtained by solving Eq. (HD and numeri¬ 
cally differentiating {E^)n^ with respect to a. Since the 
perturbation equations assume that the dark matter is 
nonrelativistic, it is best to start the numerical integra¬ 
tion of Eq. (IA18I) in the freeze-in scenario at the lowest 
temperature for which the mode lies outside the horizon. 

Finally, since Eq. (lAISdl) is not affected by the addition 
of dark matter annihilations, the superhorizon evolution 
of 0^ is given by 


^ (B3c) 

for both freeze-in and freeze-out scenarios M- When 
Eqs. (IBIIl . (|B2p . and (IB3I1 are used to set the initial con¬ 
ditions for the perturbations, the modes that remain out¬ 
side the cosmological horizon during the EMDE evolve 
to standard adiabatic initial conditions during reheating: 

—>■ 2$ = 20/9<i>o, and —>■ 3/45^- When these modes 
enter the horizon, they follow the standard evolution for 
density perturbations in a radiation-dominated universe: 
6r oscillates, while 6^ grows logarithmically. 
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